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FOREWORD 


This  collection  contains  articles 
written  by  staff  members  of  the  depart¬ 
ment  of  radio  engineering  fundamentals 
and  antenna  feed  installations  of  the 
Radio  Engineering  Division,  Moscow  Power 
Institute,  in  1971. 

As  to  subject  matter  orientation,  the 
collection  is  divided  into  four  parts. 

Part  1  includes  studies  on  problems  of  the 
electrodynamics  of  inhomogeneous  media. 
Presented  in  Part  2  are  articles  in  which 
problems  of  the  excitation  and  diffraction 
of  electromagnetic  waves  are  examined. 

Part  3  contains  investigations  on  the  re¬ 
cording  of  electromagnetic  radiation  in  the 
microwave  band.  Part  4  considers  problems 
of  the  theory  and  application  of  dielectric 
waveguides . 

Because  of  the  subject  matter  orienta¬ 
tion  of  the  collection,  the  references  cited 
are  general.  In  this  collection,  a  three- 
digit  arrangement  is  used  in  citing  articles. 
Thus,  in  referring  to  article  7,  the  number 
107  is  given,  and  In  citing  article  23,  the 
number  123. 

Much  of  the  work  Involved  in  collecting 
the  articles  and  editing  the  collection  was 
done  by  candidate  of  technical  sciences  V.  V. 
Shtykov.  Considerable  work  in  the  technical 
preparation  of  the  collection  was  performed 
by  degree-seeking  student  I.  T.  Krylova. 
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PART  1 

ELECTRODYNAMICS  OP  INHOMOGENEOUS  MEDIA 


SOLVING  PULSED  PROBLEMS  OP  ELECTRODYNAMICS  IN  INHOMOGENEOUS 
MEDIA  BY  THE  METHOD  OF  FINITE-DIFFERENCE  NETS 

Ye.  A.  Filatova 

Modeling  of  Initial  Equations  in 
the  Case  of  Magnetic  Type  Waves 

An  examination  is  made  of  modeling  with  finite-difference 
electrical  nets  of  two-dimensional  Maxwell  operator  equations 
in  spherical  coordinates  originating  in  cases  when  the  con¬ 
ductivity  c(oj)  and  the  dielectric  constant  e  (w)  of  a  plasma 

CL 

medium  vary  with  respect  to  radius  and  the  angular  coordinate  0. 
Schemes  of  net  models  are  presented  in  an  investigation  of 
waves  of  the  magnetic  or  electrical  type. 

When  problems  of  the  propagation  of  pulsed  signals  in 
dispersive  media  are  investigated,  the  initial  Maxwell  equations 
for  instantaneous  values  are  best  represented  in  the  form  of 
operator  equations.  In  the  case  under  study,  this  is  a 
system  of  scalar  operator  equations.  Here  the  system  of 
scalar  operator  equations  can  be  written  in  the  following 
form: 


fjstr  h  J  •  W. 

£  b  C  '•  />/».*•  ( P ). 


(l) 


where  £,<p  ).  /£<P  ).  H,{p  )  are  the  component  images  of  the 

vectors  of  electromagnetic  field 
intensity  after  Laplace 
m is  the  complex  dielectric  constant 
of  a  medium  in  operator  form 

The  resulting  operator  equations  are  valid,  assuming  zero 
initial  conditions. 

The  method  of  modeling  operator  equations  consists  of 
writing  the  initial -equations  in  the  form  of  equations  in 
finite  differences  for  a  constant  quantization  step  between 
the  nodes  of  a  two-dimensional  net  mn  and  mn+1  or  mn-1  and 
mn,  respectively,  for  the  step  A9  between  the  nodes  and  m+ln 
or  m-ln  and  mn.  By  introducing  the  constant  modeling  coeffi¬ 
cients  k1,  k2,  and  the  following  notation: 


,  J  (p) 

A  r 

&  tianMMUkLBlm  -  tr  tp\ 

a  9  ■*»•!**«/*, 


we  can  write  Eq.  (1)  in  finite  differences  in  the  form  of  the 
following  system: 


CL..(P)’*L.(p)'  PK&r^  sinOm  (P)t 

^L(P)’^L,  (p)  «  Pfittr  %  On 9m  J_.M  (P)t 

^  (P)  ~  am.«JP)  *  (p)* 

m*?PZ~(P)&*ZLL(p). 


(3) 


The  system  of  equations  (1)  can  be  reduced  to  operator  equations 
(8)  of  the  voltage  U(p)  balance  and  of  the  current  J(p)  balance 
in  the  node  of  the  two-dimensional  net.  The  two-port  networks 
of  the  model  can  be  calculated  with  the  following  formulas: 


(P)  mPMt  , 


Z«.,M  (P)  "  fl/vrf.  <>„  , 

Lw  s-,- 


The  model  can  be  obtained  if  from  the  resulting  equations  (4), 
(5),  (6)  the  target  is  synthesized  for  the  known  law  e  (p). 

a 

For  the  case  of  a  plasma  medium,  the  dielectric  constant 
e„(p)  and  the  conductivity  a(p)  in  operator  form  can  be  written 
as  follows: 

•  ('  -  >.<>">'  ,  ( 7 ) 

10-*^.  (8) 

where  ng  is  the  electron  concentration  in  the  plasma 

v  is  the  effective  rate  of  electron  collisions  in  the 
plasma,  1/s 

By  the  synthesis  of  the  two-port  networks,  the  model  dia¬ 
gram  as  shown  in  Fig.  1  was  obtained. 

The  model  elements  can  be  calculated  with  the  formulas: 


i  K.  ri  A  6*  ,  n 

«»  "A  Ar  » 

(9) 

(10) 

• 

(11) 

3 

i  .  *1* iV— — . 

""  **  V** 

D  m  «V  Vmm  tin  8  m 
''■»  /f.  ntri^et\F  * 


(12) 

(13) 


Modeling  Initial  Equations  in  the  Case  of  Electric  Type  Waves 


Operator  Maxwell  equations  for  this  case  can  be  written 
in  the  following  form: 

fifcl  isiHvlf&ttnO. 7-  PK(P)E*(P), 
l  &  [rHf(p)3  •  -p?Jp)E,(p), 

r  &.CrE0(ptf- 1  H^mmPf^<P)-  d4) 

Modeling  of  system  (14)  is  analogous  to  modeling  of  Eqs.  (1). 
The  model  diagram  is  given  in  Fig.  2.  The  model  elements  for 
the  case  of  a  plasma  medium  can  be  calculated  with  the  follow¬ 
ing  formulas: 


(15) 

Lmur,m  £.  ~tr~  st-nemt9r  ’ 

(16) 

’  _  K.tr 

W.  Kt€. 

(17) 

jQ  .  _  £a  X*  i,19  t0  F  4®  SiflQmtf  f 

■  -*M*  pr  Vm4M„ 

(18) 

L~nt%'U£"%*r.SLne”> 

(19) 

.  • 

(20) 

(21) 
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In  this  diagram  the  components  of  the  electromagnetic  field 
can  be  determined  from  the  following  relationships: 


(22) 

",  &&$***■  ■ 

(23) 

«.  £~7s  m-  *  J—.<~  (pi. 

(24) 

where  k^,  k^  are  constant  coefficients  of  modeling. 

The  solution  of  the  practical  problems  can  be  reduced  to 
setting  up  the  diagram  of  the  models  according  to  Pig.  1,  when 
the  required  boundary  conditions  are  satisfied  [1]. 

The  calculations  and  the  experiments  showed  that  it  is 
entirely  possible  to  exactly  synthesize  the  modeling  diagrams 
in  the  frequency  band  exceeding  the  band  used  in  actual  radio 
systems.  The  advantages  of  the  modeling  method  consist  of 
the  fact  that  in  this  case  use 'is  made  only  of  quantization 
along  the  coordinates  along  which  the  parameters  of  the  medium 
vary.  In  contrast  to  numerical  methods,  this  method  does  not 
require  quantization  with  respect  to  frequency  and  to  time. 

The  model  is  best  constructed  for  a  frequency  distinct  from 
the  working  frequency  of  the  pulse  system.  All  elements  of  the 
diagram  and  the  pulse  signals  tested  are  recalculated  to  the 
model  frequency  according  to  the  principle  of  electrodynamic 
similitude.  The  excitation  condition  is  attained  by  connect¬ 
ing  a  radio  pulse  generator  to  the  corresponding  node.  When 
the  shape  or  duration  of  the  radio  pulse  is  varied,  the  model 
elements  remain  unchanged.  A  pulse  signal  is  observed  in  the 
electrical  circuit,  at  different  points  in  it  corresponding  to 
certain  spatial  points  of  the  electrodynamic  problem. 
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SOLUTION  OF  ONE-DIMENSIONAL  PROBLEMS  OF  LONG-RANGE  ELECTRODYNAMIC 
PROBLEMS  BY  METHOD  OF  MODELING  ON  MODELS  WITH  SHORT  ELECTRICAL 
LENGTH 

I.  I.  Chugunov  and  L.  P.  Abasheyeva 

When  a  plane  electromagnetic  wave  is  incident  at  an  angle 
0  on  a  plane  inhomogeneous  layer,  the  complex  dielectric  con¬ 
stant  being  and  varying  along  the  coordinate  z,  de¬ 

pending  on  the  polarization,  the  process  can  be  described  with 
either  the  equation 

(i) 

or  the  equation 

cfa  *  6t)  %( *  m°>  ( 2 } 


where 

j(i)  •  wfaC'CEd)  -  stnVJ. 

Using  the  discrete  partitioning  of  the  layer  into  m  equal 
parts  Az,  within  whose  limits  the  variation  in  the  parameters 
of  the  medium  is  insignificant,  and  by  applying  the  principle 
of  electrodynamic  similitude,  let  us  model  Eqs.  (1)  and  (2) 
[2]. 
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By  modeling  the  equations,  we  get  the  circuit  diagrams 
shown  in  Pig.  la  and  lb. 


For  Eq.  (1) 
given  the  values 

Hence  we  have 


Fig.  la 
f (z)<0. 


,  given  the  values 
In  these  diagrams 

.  .-v  •  Ip- — s 


f(z)>Q,  and  Fig.  lb, 
according  to  [1] 


•  it**** , 

‘ 


or 


•  (5) 


For  the  case  when  ^t%g- ,  the  inductance  Lm+0  ^ 
can  be  replaced  with  the  inductance  Cm+Q  ^  determined  with 
the  formula 

c  m  K, _ 

KvM**£.(£m,i,-Un*'f)  (6) 

In  these  circuits  the  voltage  at  the  nodes  determine  the 
magnitude  of  the  magnetic  field  intensity  at  the  corresponding 
points  of  the  inhomogeneous  layer.  The  value  of  the  current 
?mtis  determines  the  electric  field  intensity.  By  modeling 
Eq.  (2),  we  get  circuits  of  the  form  shown  in  Fig.  lb,  for 
f(z)>0,  and  Fig.  Id,  for  the  values  f(z)<0. 

As  can  be  shown  in  Fig.  lc  and  Id,  in  place  of  the  con¬ 
stant  capacitance  Cm  in  the  parallel  branch  of  the  circuit 
the  parallel  connection  of  capacitance  Cm  and  resistance  Rm, 
determined  with  the  following  formulas,  now  appear: 


8 


m  At  • 

^  * 


RmtV*7<iA*  S  "'«*  ' 


*  x/  **  ( 10 ; 

For  the  negative  values  em+Q  the  inductance  Lm+Q  ^  Is  re¬ 
placed  with  the  capacitance 


frntinW 


then  the  capacitance  Cm  must  be  replaced  with  the  inductance 


In  this  diagram,  the  voltage  at  the  nodes  Um  determine 
the  value  of  the  component  of  the  magnetic  field  intensity, 
and  the  currents  Z*+  ,  the  value  of  the  electric  field  inten 
sity . 


In  Eqs .  (3)  through  (12),  and  k2  are  constant  coeffi¬ 
cients  selected  in  accordance  with  the  principle  of  electro¬ 
dynamic  similitude. 


The  design  of  the  universal  model  (one-dimensional)  makes 
it  possible  to  rapidly  select  different  diagrams  shown  in 
Pig.  1,  and  to  match  the  model  input  with  the  generator  output. 
Additionally,  there  is  the  option  of  changing  the  model  *** 
according  to  the  necessary  boundary  conditions. 

When  electromagnetic  waves  pass  through  an  inhomogeneous 
plasma  layer,  the  attenuation  constant  a  is  the  main  parameter. 

But  when  the  electrical  length  of  an  inhomogeneous  plasma  layer 
exceeds  the  electrical  length  of  the  model  or  when  the  inhomo¬ 
geneous  plasma  layer  has  large  attenuations,  U cannot  always 
be  measured.  In  these  cases,  the  problem  must  be  solved  by 
stages.  To  do  this,  the  inhomogeneous  plasma  layer  is  divided 
into  parts. 

A  block  diagram  for  determining  attenuation  by  parts  is  shown 

in  Fig.  2a.  The  letters  A  and  D  designate  the  start  and  end  of 

c. 

the  model,  respectively;  the  letters  B  and  9,  the  points  dividing 
the  inhomogeneous  layer  into  three  parts:  I,  II,  III.  When  the 
problem  was  formulated.  It  was  assumed  that  the  plane  wave  is 
incident  from  free  space  at  an  Inhomogeneous  layer  and  exits 
into  unbounded  free  space. 

As  can  be  seen  in  Pig.  2b,  the  solution  is  carried  out  to 
completion,  that  is,  first  part  III  of  the  layer  is  considered. 

The  of  the  section  of  the  model  CD  Is  determined  by  con¬ 

necting  the  generator  to  the  model  at  point  C  across  resistance 
R,  serving  for  determining  the  complex  amplitude  of  the  current, 
and  the  attenuation  constant  for  layer  III. 


By  replacing  the  section  of  the  model  CD  with  the  resistance 
*4j.*  obtained,  we  determine  ot^j  and  and  so  on.  If  the  number 

of  partitionings  is  m,  total  attenuation  is  obtained  by  summing 
the  attenuation  constants  am. 


10 


When  the  signal  passes  through  the  inhomogeneous  layer, 
the  phase  of  the  arriving  signal  with  respect  to  the  input  is 
measured  in  the  usual  way:  by  connecting  a  phasemeter  to  the 
model  input  and  output.  Shown  in  Fig.  2c  is  a  phasemeter  con¬ 
nection  diagram.  But  if  the  signal  does  not  pass  through  the 
layer,  the  above-described  method  of  dividing  the  model  into 
parts  is  used.  Then  the  total  phase  bAD  will  be 

The  phasemeter  connection  diagram  for  connecting  the  phase¬ 
meter  to  each  individual  model  part  is  the  same  as  in  Fig.  2c. 
Measuring  the  phase  by  parts  yielded  good  agreement  with  the 
result  of  measuring  the  phase  for  the  entire  model.  The  error 
of  measurement  was  +10.  Actual  working  with  the  universal 
model  showed  that,  by  successively  replacing  the  model  parts 
with  their  >  because  of  rapid  commutation  of  the  model 
elements  the  problem  can  be  solved  without  significant  degrada¬ 
tion  of  solution  accuracy. 
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REFLECTION  OF  A  PLANE  TM  WAVE  FROM  A  HOMOGENEOUS  PLASMA 
HALF-SPACE 


I.  G.  Yakushkin  and  V.  A.  Permyakov 

1.  As  is  known,  when  a  TM  wave  is  reflected  from  an  in¬ 
homogeneous  plane  stratified  plasma  with  the  law  of  variation 
of  the  dielectric  constant  >  in  the  neigh¬ 

borhood  of  zero  e,  resonance  absorption  of  the  field  occurs. 

In  this  case,  the.  modulus  of  the  reflection  coefficient  is 
different  from  unity  even  when  there  are  infinitesimal  heat 
losses  in  the  layer.  Although  the  phenomenon  of  resonance 
absorption  has  been  investigated  in  a  large  number  of  publica¬ 
tions  [3-6],  they  lack  detailed  numerical  results  and  approxi¬ 
mate  expressions  suitable  for  determining  the  reflection  coef¬ 
ficient  when  there  are  arbitrary  gradients  of  the  dielectric 
constant.  Below  we  analyze  the  function  of  the  reflection  coef¬ 
ficient  of  a  TM  wave  reflected  from  a  plasma  half-space,  with 
the  linear  law  £(»)•-«■,  D-jf-  when  there  are  arbitrary  gradients 
of  e . 


The  reflection  coefficient  for  the  magnetic  field  is  defined 
by  the  expression 


»  cot9~ 
cost 


Sf 


(1) 


where  Z  is  the  normalized  (divided  by  free-space  impedance) 
impedance  of  the  layer,  equal  to 


(2) 


,  .  [/'(*) 

Z  *  **•  K£(i>Ull) 

and  U(z)  satisfies  the  equation 

The  angle  of  incidence  of  the  wave  9  is  measured  from  the 
normal  to  the  interface. 

2.  When  there  are  arbitrary  gradients  of  e,  the  solution 
of  Eq.  (3)  is  sought  for  by  the  power  series  method  at  the 
point  e*0  and  has  the  form  [4,  9] 

uw  •(/«  +  t/aah  (4) 

where 

ujn-Ia.i:  a,-a,.o. 

The  unknown  function  f(A)  appearing  in  Eq.  (4)  is  determined 
on  the  condition  that  the  solution  of  U(t)  decreases  as  t-M-°° 

Different  formulas  are  presented  in  [9]  for  determin¬ 
ing  the  functions  f(A)  suitable  for  both  analytic  and  machine 
calculations.  According  to  [ 9 ] ,  when  |A|<1,  the  function  f(X) 
is  equal,  with  at  least  an  accuracy  of  5  percent,  to  its  value 
when  9=0:  /(*)*/(«&*  0,8M. 

Using  this  value  of  the  function  f(A)  and  retaining  the 

p 

terms  of  the  power  series  (4)  to  the  order  (k/a)  inclusively, 
we  get  the  following  asymptotic  formula  for  skin  impedance  that 
is  suitable  for  large  gradients  of  e: 

?  .  *19 t]$  Sin1* 
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Beginning  with  Eq.  (5),  we  can  show  that  when  a>>l,  the  modulus 
of  the  reflection  coefficient  takes  on  the  minimum  value  of 
|R|rain%1  to  2.3(k/a)2//^  when  e^arccos  ImZ.  With  increase  in  a 
Rmin  tends  to  unity  and  is  observed  at  grazing  angles  of  inci¬ 
dence  . 


3.  When  a<<k  and  for  not  very  shallow  incidence,  the  phase 
integral  method  [7]  can  be  applied  for  an  approximate  solution 
of  Eq.  (3).  The  phase  Integral  method  solution  of  Eq.  (3)  can 
be  constructed  in  both  the  region  of  large  positive  z,  that  is, 
where  e<0,  as  well  as  in  the  region  of  negative  z.  After  the 
substitutions  ,  Eq.  (3)  takes  on  the*  form 

(6) 

where 

A  tin*6 

For  large  in  Eq.  (6)  the  term  ( 3/^ ) can  be  neglected. 
So  the  phase  integral  method  solution  of  Eq.  (6)  for  large  nega¬ 
tive  £,  selected  on  the  conditions  of  decreasing  in  the  domain 
where  e<0,  takes  on  the  form 

-frr ^  (7) 

For  large  positive  C,  based  on  the  radiation  conditions,  we  can 
write 

f .  (np(-i J fc-A  dg  +R9xp(if\/^Td$)]  (8) 

The  coefficient  A,  in  the  case  if  we  can  neglect  reflection 
from  the  interface  of  media,  can  be  determined  simply  on  the 
condition  of  interlinking  with  the  incident  wave.  The  reflection 
coefficient  of  the  wave  R  from  an  inhomogeneous  half-space  is  the 
principal  quantity  to  be  determined.  For  large  values  of  A, 
generally  the  energy  does  not  penetrate  to  the  point  £=0.  In 
this  case,  Eq.  (6)  can  be  reduced  to  the  Airy  differential  equa¬ 
tion  and  it  turns  out  that  R-l.  The  situation  is  more  involved 
when  A^l,  since  in  this  case  the  reflection  is  determined  by  the 
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function  g-A-£g*4  .  This  function  has  already  three  inflection 
points,  two  of  which  are  located  at  complex  values  of  £. 


To  take  this  reflection  into  account,  by  using  the  phase 
Integral  method  it  is  important  to  be  able  to  consider  the 
phase  integral  method  approximation  to  be  applicable  even  in 
the  neighborhood  of  a  pole.  For  this  purpose,  the  function 
p(£)  used  in  constructing  the  phase  integral  method  solution 
[71,  Eqs.  (6)  must  be  selected  so  that  the  error  of  the  phase 
integral  method  solution  remains  finite  when  £=0.  The  selection 
[6]  of 


/xe)  \g  -A-f-* 


(9) 


leads  to  this  result. 


j.n  the  following  treatment,  the  notation  below  will  be  used 

CS,*3-"P(-cfip(s)Js).  *  .  do) 

t 

To  apply  the  phase  integral  method  in  determining  the  coupling 
coefficient  associating  the  asymptotic  representations  of  the 
solution  for  positive  and  negative  £,  we  must  know  the  structure 
of  the  Stokes  lines,  that  is,  the  lines  where 


(£  is  the  root  of  the  equation 
p(£)=0). 


The  pattern  of  the  Stokes  lines  for  the  function  p(5)  assigned 
by  Eq.  (9)  is  shown  in  Fig.  1. 

Jugate  Stokes  lines  at  which 


The  dashed  line  shows  the  con- 


For  large  negative  £>  the  solution  must  be  sought  in  the 
form  of  (7),  that  is,  in  region  IV  the  solution  is  of  the  form 
F*CU,52). 
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Pig.  1 


In  region  III  this  solution  is  growing  and  therefore  in  the 
transition  to  region  II,  in  accordance  with  the  conventional 
rules  for  transition  through  conjugate  Stokes  lines  [7,8],  we 
get 


At  Stokes  line  I  we  have 

F-C{<1  or S.SJ  ♦  Its.  ♦  iSjXlJj}. 

By  selecting  C  so  that  the  solution  is  brought  to  the  form 
(8),  that  is,  ,  where 

exp  ♦£ /  vjp-A  dg}t 

*  ^  • 

we  get 

ft  •  l9r*{i*Cgatgf]aJ  (ID 


To  derive  the  formula  for  the  reflection  coefficient,  we 
used  the  complex  inflection  point  lying  in  the  upper  half-space. 
Use  of  the  lower  inflection  point  does  not  lead  to,  and  must 
not  lead  to  the  same  result,  since  Eq.  (4),  at  point  C=0,  has 
a  branch  point  of  infinite  order,  therefore  by  bypassing  this 
point  in  the  lower  half-space,  we  arrive  at  another  sheet  of  a 
manifold  plane  on  which  the  solution  is  defined. 

The  cofactor  appearing  in  Eq .  (11)  ahead  of  the  brackets 
for  the  modulo  is  equal  to  unity  and  describes  the  phase  shift 
of  the  reflected  wave.  The  cofactor  in  the  brackets  describes 
the  resonance  absorption  of  the  TM  wave. 

In  the  limiting  cases,  the  solution  obtained  coincides  with 
known  results. 

As  can  be  seen  from  (11),  R+i  as  A-*-».  When  A=0,  the  same 
value  R=i  is  obtained  from  the  exact  solution  and  from  Eq.  (11). 


.  .  -m 


4.  Pigs.  2  and  3  present  the  results  of  computing  the  re¬ 
flection  coefficient  by  the  power  series  method  (4)  on  a  com¬ 
puter  (solid  line),  and  their  comparison  with  the  approximation 
of  a  strongly  inhomogeneous  plasma  (5)  (crosses)  with  the  appro¬ 
ximation  of  a  weakly  inhomogeneous  plasma  (dashed  line)  (11). 

From  a  comparison,  there  must  be  close  agreement  of  the  approxi¬ 
mation  of  the  weakly  inhomogeneous  plasma  with  the  rigorous  cal¬ 
culation  when  a<_0.1k  in  the  entire  range  of  angles  of  incidence 
for  plane  waves.  Prom  the  graphs  it  follows  that  j  R |  has  a  min¬ 
imum  equal  to  ^0.7  when  X^0.5j  and  beginning  with  \^2.25,  |R|  is 
in  practice  equal  to  unity.  For  large  gradients  of  e  (l>a/k> 
>0.1),  the  phase  integral  method  is  suitable  in  the  range  of 
angles  close  to  zero  (9<40  to  60°).  The  inaccuracy  of  approxima¬ 
tion  (11)  at  grazing  angles  of  incidence  in  this  region  of 
gradients  of  e  is  explained  by  the  fact  that  in  its  derivation 
the  layer  is  assumed  unbounded,  while  the  numerical  results  were 
obtained  when  the  wave  was  incident  on  a  half-space. 


1 


The  approximation  of  a  strongly  inhomogeneous  plasma  (5) 
proves  to  be  suitable  when  a^k.  It  bears  noting  that  when 
a^(l  to  0.5)k  in  the  region  of  grazing  angles  of  incidence, 
this  approximation  agrees  closely  with  the  rigorous  calculation. 

Thus,  in  the  aggregate  the  approximate  solutions  for  the 
weakly  inhomogeneous  and  strongly  inhomogeneous  plasma  make  it 
possible  to  describe  the  behavior  of  the  reflection  coefficient 
of  the  TM  wave  for  arbitrary  gradients  of  e. 
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EFFECT  OF  LINEAR  TRANSITION  LAYER  ON  THE  ATTENUATION  OF  WAVES 
IN  A  PLANE  STRATIFIED  WAVEGUIDE 


V.  K.  Polishchuk 


1.  Derivation  of  Dispersion  Equations 


Let  us  examine  a  plane-stratified  waveguide,  with  relative 
dielectric  constant  er,  which  is  a  continuous  even  function  of 
the  coordinate  x  (Fig.  1)  and  which,  when  x>0,  can  be  expressed 
by  the  formulas 


£„(*) 


r o*k<(, 

**d, 


If  a  slit  in  a  screen  symmetrically  positioned  and  symme¬ 
trically  excited  with  respect  to  the  waveguide  axis  is  excited 
with  a  source,  only  even  TM  waves  will  exist  in  the  waveguide  [10]. 


There  is  a  single  component  of  the  magnetic  field  H  in  the 
TM  wave  propagating  along  the  z  axis.  Considering  that  the  de¬ 
pendence  of  H  on  coordinates  is  of  the  form 

(we  assume  that  the  field  does  not  depend  on  the  coordinate  y) 


Fig.  1 


h-*£h'h«X-i,‘)h  ■».  a, 

where  kQ  Is  the  free-space  wave  number. 

The  component  of  the  electric  field  E  that  is  tangent  to 
the  layer  boundaries  can  be  determined  in  terms  of  H(x)  with  the 
formulas 


E(x)  »  (iu>£,£r)'tH  (x). 

where  e Q  is  the  absolute  dielectric  constant  of  a  vacuum.  Let 
us  find  the  expressions  for  the  functions  H(x)  and  E(x)  in  dif 
ferent  regions  of  the  waveguide;  in  view  of  symmetry,  we  will 
consider  the  field  only  when  x>0. 

Region  I.  0<x<£ 

Eq.  (1)  takes  on  the  form:  H"+q^H=0.  Solving  it,  for  the 
even  wave  types  we  get 

H(x)  •  cos  $x, 

?  ♦  h‘.  . 
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The  expression  for  H(x)  is  determined,  strictly  speaking,  with 
an  accuracy  to  an  arbitrary  constant  cofactor.  But  since  the 
amplitude  of  the  waves  in  the  waveguide  is  not  of  interest  to 
us,  we  assume  that  in  region  1  it  is  equal  to  unity,  for  H(X). 

Region  3.  x>£+d 

Eq.  (1)  takes  on  the  following  form:  H"-p^H=0.  With  allow 
ance  for  the  weakening  of  the  field  as  x-*°°,  we  find 

£(*)• p(tiae,tS K*’*  ( . } 

p'-h'.ptt.. 

Region  2  [sic].  £<x<J!.+d 

In  region  2  we  have  the  equation 

H“*  jfciH'+fafa (x,-x)-h*)H  »o  ( 4 ) 

We  will  solve  it  using  series  in  powers  of  the  difference 
xQ-x  and,  considering  the  width  of  region  2  to  be  small 
(d/2.<<l),  in  the  result  we  will  limit  ourselves  to  a  finite 
number  of  terms. 


It  is  known  [11]  that  the  general  solution  of  Eq.  (4)  is 


of  the  form 


H(x)  •  ctu(x)*ctv(x). 


where  u  and  v  are  two  linearly  independent  solutions  of  Eq.  (4) 
represented  in  the  form  of  the  expansions 

U  V,  suln K.(X.*X) 

whose  coefficients,  determined  from  the  recursion  formulas, 
are  equal  to,  respectively: 

4  •  £ »  4>  *  V  A*  •  ■  & . 
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Using  these  expressions,  by  limiting  ourselves  to  terms  not 
higher  than  the  first  order  of  smallness  in  Xq-x  and  d,  we 
get 

*At  E(x)  •  •(itot'O)  (Ah 1  In  k,(x,-x)  *£)  ( 5 ) 

By  here  equating  the  corresponding  Eqs.  (2),  (3),  and  (5)  at 
the  boundaries  of  regions  1  through  3,  we  get  the  following 
system  of  dispersion  equations: 

9  if  h  *&)* 

p*+  q,\m  *}(£.+  O  (6) 

2.  Analysis  of  Dispersion  Equations 


By  neglecting  waves  of  all  types  except  for  the  principal 
wave,  we  analyze  system  (6)  for  the  symmetrical  transition  layer 


:+=e_=e. 


If  we  convert  to  the  quantities  , 

and  ,  and  if  we  use  Eqs.  (2)  and  (3),  we  can  get  the  follow¬ 

ing  system  of  eauations: 

9  *xtax  • 

]  (7) 

cr  the  single  dispersion  equation  for  the  longitudinal  wave 


number 


To  find  the  solutions  of  the  system  (7)  or  Eq.  (8),  we  assume 
that  they  are  known  to  use  when  6  =  0  and  .  We  denote  the 

quantities  a,  t,  and  p  corresponding  to  this  case  with  Pq,  tq, 
and  Pq,  that  is,  under  our  assumption  the  following  relationships 
obtain: 


.«.*T 

3Li£j7**lL 

Vir •© 


Suppose  now  £  =  ,  and  6^0.  Let  us  represent  the  longi¬ 

tudinal  wave  number  PQ  in  the  form  papg+p^  and  substitute  it  in 
Eq.  (8).  By  expanding  the  left  and  right  parts  of  (8)  in  series 
in  terms  of  the  small  parameters  6,  e^,  and  and  referring  to 


Eq.  (9),  we  get 


J>  •  f>J<  -  t(*S 

with  an  accuracy  to  the  infinitesimals  of  the  second  order,  where 


We  see  that  to  the  first  approximation  of  the  method  of  per¬ 
turbations,  the  corrections  to  the  constant  of  propagation  due  to 
the  transition  layer  and  to  the  imaginary  part  of  the  dielectric 
constant  are  additive  and  therefore  can  be  examined  independently 
of  each  other.  In  addition,  in  the  case  of  a  symmetric  transi¬ 
tion  layer  these  corrections  are  purely  imaginary  for  propagating 
wave  types,  that  is,  their  phase  velocities  remain  unchanged. 


Presented  in  the  study  [12]  are  the  solutions  of  the  "un¬ 
perturbed"  equation  (10)  for  different  »  .  When  investigating 
the  behavior  of  the  correction  coefficients  a  and  3,  it  is  con¬ 
venient  for  us  to  use  Eq.  (10)  as  well  as  system  (9).  And  as  the 
independent  variable  we  select  not  the  frequency  « ,  but  the 
transverse  wave  number  tq. 


Let  us  find  the  critical  values  of  ar  determining  the  reso¬ 
nance  in  the  cross-section  of  the  waveguide  at  which  Pq=0. 


By  substituting 
from  which  we  find 


p  =0  in  Eq.  (10),  we  get  the  equation 
0 

1*0 *Ot 


I  0,  m»0 

m*0. 


By  differentiating  Eq.  (10)  at  the  point  «•  *f,  we  can  show  that 
when  the  longitudinal  wave  number  pQ  is  imaginary,  and 

when  *>s*»  it  is  real,  that  is,  the  critical  value  of  for  the 
fundamental  wave  is  ^  ,  and  at  *.-0.  If  the  frequency  «  takes 
on  the  values  within  the  limits  *<*<?  ,  all  waves  of  the  higher 
types  will  be  attenuating. 
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When  %<x<  ,  the  longitudinal  wave  number  Tq  takes  on 

values  in  the  interval  &  . 

The  solutions  of  Eqs.  (9)  and  (10)  in  the  interval 
(for  the  fundamental  wave)  are  shown  in  Pig.  2  in  the  form  of 
the  functions  <»,(*,)  •  ae.fr;)  and  .  Also  presented  there  are 

the  plots  for  the  coefficients  “(Tq)  and  8(tq)  determining  the 
behavior  of  the  imaginary  part  of  the  constant  of  propagation 
P,  that  is,  the  decrement  of  attenuation.  From  Pig.  2  we  see 
that  for  the  fundamental  wave  the  curve  8(tq)  trends  everywhere 
higher  than  the  curve  oi(Tq),  such  that  when  there  are  identical 
values  of  6  and  e^,  the  resonance  absorption  will  always  be  less 
than  the  heat  losses.  In  addition,  when  r^-*-Tr/2,  ae— (at  higher 
frequencies)  the  function  S(tq)  tends  to  a  finite  value,  and 
a(t3),  to  zero.  So  by  increasing  the  frequency  (or,  which  amounts 
to  the  same  thing,  the  electrical  dimension  of  the  waveguide),  the 
resonance  absorption  can  be  made  small  enough. 

In  the  other  limiting  case,  when  the  frequency  is  near- 
critical  in  value,  we  can  find 

3-JT  4Z3T  i 

f  T/' 

In  this  case  we  observe  the  strong  dependence  of  8  on  fre¬ 
quency  and,  in  addition,  8-”°  (t-*t  ) .  So  at  near-critical  fre¬ 

quencies,  Eq.  (11)  may  prove  to  be  invalid  for  8  obtained  by  the 
method  of  perturbations. 

Calculations  show  that  in  estlmational  calculations  of  the 
attenuation  of  radio  waves  in  an  inhomogeneous  waveguide  the 
losses  due  to  the  transition  region  can  be  neglected. 
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CHARACTERISTICS  OF  RADIATION  OF  APERTURE  ANTENNAS  IN  A 
PLANE-STRATIFIED  MEDIUM 


A.  ?.  Anyutin,  V.  A.  Permyakov,  and  V.  K.  Polishchuk 

Presented  in  this  study  are  working  formulas  and  the  result 
of  calculation  of  radiation  patterns  of  a  (circular)  waveguide 
located  in  an  infinite  metal  plane  under  a  dielectric  layer  and 
a  plasma  layer.  In  deriving  the  formulas  for  the  radiation 
patterns  it  is  assumed  that:  1)  the  field  distribution  in  the 
open  end  of  the  rectangular  (circular)  waveguide  coincides  with 
the  field  distribution  of  the  wave  HQ1  (H.^)  propagating  in  the 
waveguide;  2)  the  reflection  coefficient  p  of  the  wave  HQ1  (H-^) 
from  the  open  end  of  the  waveguide  Is  given;  and  3)  the  complex 
dielectric  constant  of  the  plasma  varies  only  along  the  normal 
to  the  surface  and  is  defined  by  the  formula  [3]: 


where  w  is  the  working  frequency,  vg  is  the  effective  number  of 
electron  collisions  in  the  plasma,  and  Ng(z)  is  the  electron  con 
centration  in  the  plasma.  These  assumptions  allow  us  to  reduce 
the  initial  problem  to  a  problem  of  the  radiation  of  a  current 
sheet  In  a  two-sheet  medium  (Fig.  1).  A  method  of  solving  this 
problem  consists  in  the  following  [13]:  by  applying  the  Fourier 
transform  to  the  wave  equation  with  respect  to  the  coordinates 
x,  y,  and  by  satisfying  the  boundary  conditions,  a  system  of 


linear  algebraic  equations  results  for  the  spectral  representa¬ 
tions  of  the  fields.  The  dependence  of  all  spectral  representa¬ 
tions  on  the  coordinate  z  is  determined  with  ordinary  differen¬ 
tial  equations  with  variable  coordinates  (in  the  case  of  an  in¬ 
homogeneous  plasma). 


Pig.  1 

By  applying  the  inverse  Fourier  transform  to  the  solution  of 
the  system  of  linear  algebraic  equations,  we  get  the  integral 
representation  of  the  fields.  To  define  the  fields  in  the  far 
zone,  the  double  integrals  representing  them  are  calculated 
with  the  method  of  the  integral  [14],  The  final  expressions 
for  the  components  of  the  electric  field  in  the  far  zone  are 
as  follows: 

a)  rectangular  waveguide 


-  ■ 


b)  circular  waveguide 

_  r«  *  Plr.  cot  a  »i  fv-p)  s  ine)cose  un  t 

In  Eqs.  (1)— (4) ,  U(z)  is  the  solution  of  the 

equation 

+  *:&,(*)  -  Sin'MlU'O 

(5) 

with  the  initial  conditions  U(4)m  I  .  W(t,)*-iK.co%0 

point  z-i2>  V(z)  is  the  solution  of  the  equation 

at  the 

int  St*  K^[£t(M)-UffteJV’0 

(6) 

with  the  initial  conditions  I,  V'(ti)—u<tej(£,)coiO  at  the 

point  z=&2> 

a  is  the  dimension  of  the  narrow  wall  of  the  waveguide 
(waveguide  radius) 

b  is  the  dimension  of  the  wide  wall  of  the  waveguide 
e-^(z)  is  the  relative  dielectric  constant  of  the  plasma 
is  the  relative  dielectric  constant  of  the  dielectric; 

ft 

kQ  is  the  wave  number  of  free  space; 

) 

are  Bessel  functions 

are  as  shown  in  Pig.  1. 
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The  solutions  of  Eqs.  (5)  and  (6)  are  sought  for  numerically 
by  the  Runge-Kutta  method  [153;  where  e^(z)  is  approximated  with 
the  function  Specifying  e^(z)  in  another  analytic 

form  poses  no  fundamental  difficulties  in  calculations. 

The  expressions  found  for  the  fields  in  the  far  zone  can  be 
used  in  calculating  the  antenna  efficiency  n: 

?  “  Px/p0 » 

where  is  the  antenna  power, 

PQ  is  the  power  fed  to  the  waveguide; 
a)  a  rectangular  waveguide 


b)  a  circular  waveguide 

These  formulas  for  calculating  the  characteristics  of 
antenna  radiation  and  antenna  efficiency  in  the  plane-strati¬ 
fied  medium  were  programmed  for  computations  on  a  computer. 

For  Aq=6 . 28 ;  *2-0;  Jlj-1;  A0=4°;  a=0.03XQ;  0=0. 5AQ;  d=0; 
shown  in  Fig.  2  are  the  radiation  patterns  of  an  antenna  in 
the  plane  of  an  electric  vector  without  allowing  for  plasma 
(dot-dashed  line)  and  for  two  laws  of  variation  in  the  dielec¬ 
tric  constant  of  the  plasma,  the  linear  (solid  line)  and  the 
quadratic  (dashed  line).  The  calculated  distributions  e(z)  are 
also  presented  in  Fig.  2.  From  the  figure  it  follows  that  the 
plasma  layer  leads  to  an  abrupt  drop  in  the  relative  radiation 
level  in  the  region  of  grazing  angles  of  observation.  We  can 
also  note  the  weak  dependence  of  the  radiation  pattern  on  the 
law  e(z)  of  the  plasma,  characteristic  of  thin,  compared  with 


APPROXIMATE  METHOD  OF  DIAGNOSTICS  OF  INHOMOGENEOUS  PLASMA 

M.  B.  Aksenov,  V.  G.  Kartashev,  A.  P.  Mayorov,  and  V.  L.  Skachkov 

Theoretical  Analysis  of  Problem 

In  investigating  plasma  parameters  there  often  appears  the 
necessity  of  determining  the  law  of  spatial  distribution  of  the 
main  plasma  parameters:  concentration  Ng  and  collision  rate 

Known  methods  of  microwave  diagnostics  of  inhomogeneous 
plasma  (for  example,  [16]),  permitting  the  determination  of  the 
spatial  distribution  of  concentration,  are  based  on  the  geometric 
optical  approximation  and  therefore  are  not  always  applicable. 
Rigorous  methods  of  diagnostics  (for  example,  [17])  are  very  com¬ 
plex  and  have  not  yet  been  fully  worked  out.  But  in  some  cases, 
when  there  is  apriori  information  on  the  kind  of  distribution  of 
plasma  concentration,  the  problem  of  diagnostics  can  be  consi¬ 
derably  simplified.  Examined  in  this  study  is  a  method  of 
diagnostics  of  a  plane-stratified  plasma  based  on  the  assumption 
that  the  plasma  concentration  varies  in  space  according  to  an  ex¬ 
ponential  law,  but  the  collision  frequency  is  constant.  The  method 
is  based  on  measuring  the  coefficient  of  reflection  from  the  plasma 
when  it  is  irradiated  with  a  plane  electromagnetic  wave. 

Let  us  examine  the  normal  incidence  of  a  plane  electromag¬ 
netic  wave  from  a  homogeneous  dielectric  with  dielectric  constant 
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£-[_  on  a  plane-stratified  plasma  with  a  concentration  that  varies 
in  the  direction  of  the  z-axis  according  to  an  exponential  law 

K  (•><>).  (1) 

The  plane  z=0  is  the  Interface  of  plasma  and  the  dielectric. 
Without  constraints  on  generality,  we  will  assume  that  the  irrad¬ 
iating  field  has  the  components  E  and  H  . 

x  y 


As  shown  in  [3],  the  equation  for  the  electric  field  inten¬ 
sity  in  this  case  can  be  represented  in  the  following  form 


jfif  *  c 


(2) 


where 


_ _ _ .•  . 


(3) 


The  general  solution  of  Eq.  (2)  with  reference  to  (1)  and  (3) 
can  be  written  in  the  form 


where 


(4) 


a 

K‘ 


is  the  dimensionless  coordinate 

✓  of? 

f  *  **T  is  the  dimensionless  constant  of  attenuation 


u)m 


is  the  plasma  frequency  when  z=0. 
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From  inspection  of  the  solution  (4)  for  small  values  of  the 
arguments  u  (?“■*■)  it  follows  that  and  therefore, 

c2=0.  Finally,  solution  (4)  can  be  written  in  the  form 

e 

From  the  boundary  conditions  when  z=0  let  us  find  an  expres 
sion  for  the  complex  reflection  coefficient  F 

Using  the  recurrence  relationship  for  the  Bessel  functions,  we 
get  a  final  expression  for  the  reflection  coefficient 

r-  'S ~ M A* (a>  (5) 


Fig.  1 

By  representing  the  3essel  functions  in  the  form  of  series, 
after  some  transformations  we  get 


By  applying  Eqs.  (5)  and  (6),  we  can  calculate  the  values 
of  the  modulus  and  the  phase  of  the  coefficient  of  reflection 
by  the  plasma  for  any  value  of  the  coefficient  NQ  and  any  value 
of  the  parameter  a.  Some  results  of  the  calculations  for  e=l 
are  presented  in  Figs.  1  and  2,  where  for  convenience  in  further 
application,  the  lines  of  the  modulus  level  and  the  phase  of  the 
reflection  coefficient  are  presented  for  the  frequencies  of  the 
probe  field  f1=iO  GHz  and  f2=30  GHz. 

Method  of  Determining  the  Plasma  Parameters  Based  on  Experi¬ 
mental  Data 

Let  us  examine  the  problem  of  the  practical  determination  o 
the  plasma  parameters  Nq,  a,  and  from  measurements  of  the 
reflection  coefficient. 


If  the  collision  rate  is  small,  so  that  the  effect  of  colli¬ 
sions  on  the  reflection  coefficient  can  be  neglected,  to  deter¬ 
mine  the  plasma  parameters  Nq  and  a,  it  is  sufficient  to  measure 
the  modulus  of  the  reflection  coefficient  Ifil  and  irj  at  two 
rates  f^  and  f9.  Then  from  the  plot  shown  in  Fig.  1,  we  can 
find  the  points  of  intersection  of  the  corresponding  curves 
w  and  IV  and  determine  the  values  of  and  a. 

When  there  are  collisions  in  the  plasma,  the  procedure  of 
interpreting  the  experimental  data  is  considerably  more  involved. 
In  this  case,  to  determine  the  three  plasma  parameters  NQ,  a, 
and  ,  we  must  measure  three  parameters  of  the  reflected 

waves.  This  can  be  either  the  modulus  of  the  reflection  coef¬ 
ficient  at  three  rates  10  ,  Ifil  ,  and  iV  ,  or  the  modulus  /V 
and  the  phase  0.^  of  the  reflection  coefficient  at  the  same 
rate  and  the  modulus  of  the  reflection  coefficient  at  the 

second  rate.  Analysis  shows  that  preference  must  be  given  to 
the  second  method,  since  the  task  of  measuring  the  coefficient 
of  reflection  from  one  part  of  the  plasma  at  three  rates  proves 
to  be  complex. 

For  the  practical  determination  of  the  plasma  parameters 
in  this  case,  we  must  calculate  and  plot  the  graphs  cf  the  lines 
of  the  level  f/Jj •  »  lf\l  (Fig.  2)  for  different  values  of  a, 

thus  constructing  an  album  of  curves,  each  page  of  which  corres¬ 
ponds  to  a  certain  a.  Further,  from  the  measured  values  of  mt, 

,  and  w  we  seek  the  page  in  the  album  where  all  three 
curves  of  the  levels  |f,lt  irtl  intersect  at  the  same  point.  The 
coordinates  of  this  point  define  the  values  of  Nn  and  ,  and 
the  page  number,  the  value  of  a. 

This  method  of  diagnostics  of  an  inhomogeneous  plasma 
proves  to  be  quite  simple,  but  it  is  suitable  only  when  the 
plasma  concentration  exhibits  an  exponential  distribution.  If 
the  distribution  of  the  concentration  differs  from  the  exponen¬ 
tial,  rigorous  methods  of  plasma  diagnostics  must  be  employed. 
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PROBLEM  OF  MEASURING  THE  PARAMETERS  OF  LOSSY  INHOMOGENEOUS 
DIELECTRICS 


V.  G.  Kartashev  and  A.  P.  Mayorov 

Measurement  of  Parameters  of  Dielectrics  in  Normal  Probing 


Let  us  examine  the  normal  incidence  of  a  plane  electromag¬ 
netic  wave  on  the  surface  of  an  inhomogeneous  plane-stratified 
dielectric  with  relative  dielectric  constant  e(z).  We  assume 
the  relative  permeability  of  the  dielectric  p  to  be  constant. 

We  take  the  direction  of  the  coordinate  axes  and  the  polariza¬ 
tion  of  the  irradiating  wave  such  that  the  boundary  of  the 
dielectric  coincides  with  the  z=0  plane,  and  the  field  has  only 
the  components  Ex  and  Hy. 


(3) 

where 

z-^r;  pW'rifc&y  (4) 

Employing  the  substitution  “■m*.  ,  we  get  the  final  equation 

V  *  C.ui*p(t)  -  4(t)1C7  •o)  ( 5 ) 

where 

If  the  dielectric  under  study  is  not  lossy  (cr=0),  then 
p(x)=l,  and  Eq.  (5)  takes  on  the  form  for  which  the  inverse 
Sturm-Llouville  problem  was  solved  in  '[18].  In  this  case,  the 
problem  of  measuring  the  dielectric  constant  e(z)  of  a  nonlossy 
dielectric  can  be  reduced  to  the  following: 


The  complex  reflection  coefficient  r(to)  Is  measured  as  a 


function  of  the  field  frequency;  then  the  characteristic  func¬ 


tion  of  Eq. 


(5)  is  determined: 


4+f(w) 
cut  V-fH 


Then,  according  to  [19],  the  spectral  function  p(X)  is  found: 

.  A*# 

X)  •  tin  tin  £  /  Jm  ff(<*l*ix)du) 

and  using  the  method  given  in  [18]  the  unknown  coefficient  <p(  r) 
of  Eq.  (5)  is  determined.  Then,  by  solving  Eq.  (6),  using  Eqs . 
(U)  and  (2)  we  can  find  the  desired  function  e(z).  Specific 
procedures  for  solving  this  problem  as  applied  to  Inhomogeneous 
transmission  lines  are  considered  in  detail  in  [20] . 


If  the  conductivity  of  the  dielectric  cannot  be  neglected, 

this  method  of  solving  the  inverse  problem  is  not  applicable. 

2 

In  Eq.  (5)  the  spectral  parameter  u>  is  not  separated  from  the 
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function  of  t;  the  solution  of  the  inverse  problem  for  equations 
of  this  kind  is  not  yet  known. 

Thus  normal  probing  of  an  inhomogeneous  dielectric  with  a 
plane  electromagnetic  wave,  using  now  available  mathematical 
methods,  makes  it  possible  to  determine  the  dielectric  parameters 
only  when  the  dielectric  is  nonlossy. 

Measuring  the  Parameters  of  Dielectrics  in  Slant  Probing 

The  slant  incidence  of  a  plane  electromagnetic  wave  on  a 
plane-stratified  medium  is  described  with  the  equation  [21] 

F *  ♦  *2  -  £,  %ina0f]F  »o  ( 7 ) 

where  e'(z)  is  the  complex  dielectric  constant  of  the  test 
medium  (z>0) 

e-j^  is  the  dielectric  constant  of  the  medium  for  [text 
is  missing] 

9^  is  the  angle  of  incidence  of  the  wave 

Eq.  (7)  can  be  converted  to  the  form 

F'+U-l(*flF  • o , 

where 

A*  (* -£,Sin*Q,) 

If  the  frequency  w  is  fixed,  the  spectral  parameter  X  and 
the  function  of  the  coordinate  q(z)  are  not  related  in  Eq.  (8) 
as  with  the  real  and  the  complex  function  q(z).  The  method  of 
solving  the  inverse  problem  for  Eq.  (8)  with  the  complex-valued 
coefficient  q(z)  was  worked  out  in  [22],  To  carry  out  this 

• 

method  we  must  determine  the  complex  reflection  coefficient  r  as 
a  function  of  the  angle  of  incidence  91  and  find  the  function 
m(X) 


(8) 

(9) 

(10) 
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m(A)* 


T4cot6i  y- 


which  acts  as  the  characteristic  function  of  Eq.  (8). 

One  feature  of  this  method  is  the  fact  that  when  the  angle 
of  incidence  9,  varies ,  the  spectral  parameter  X  in  accordance 

i  2  2 

with  Eq.  (9)  can  vary  over  the  limits  from  0  to  w  /c  .  Thus, 
the  function  m(A)  proves  to  be  specified  not  for  the  entire 
semiaxis  0<A<°°,  but  only  for  the  interval  (0,oj  /c  ).  But  by 
suitably  selecting  the  field  frequency  oj  and  using  the  asymp¬ 
totic  relationships  for  m(A),'the  error  caused  by  the  incomplete 
specification  of  the  function  m(A)  can  be  reduced  to  zero. 

A  second  feature  of  the  slant  probing  method  is  caused  by 
the  displacement  of  the  beam  of  electromagnetic  waves  during 
reflection,  leading  to  a  reduction  in  the  measured  reflection 
coefficient  and  to  loss  of  information  about  the  deep-lying 
layers  of  the  dielectric.  To  reduce  the  influence  of  this 
effect,  it  is  best  to  irradiate  the  test  object  with  a  broad 
beam  whose  transverse  dimensions  are  much  larger  than  the 
assumed  beam  displacement. 
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PART  TWO 

EXCITATION  AND  DIFFRACTION  OF  ELECTROMAGNETIC  WAVES 
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SOME  SOLUTIONS  OF  MAXWELL'S  EQUATIONS 
0.  T.  Markov 

Presented  at  the  radio  engineering  section  of  the  scientific- 
engineering  conference  of  the  Moscow  Power  Institute  in 
April  1S70 

1.  Suppose  that  in  a  spherical  region  with  radius  a,  with 
parameters  of  the  medium  a=0,  eQ,  uQ  (vacuum),  there  is  given 
the  distribution  of  the  volume  density  of  electric  currents 
/*  varying  with  time  according  to  the  law  .  Imagine  an 
electromagnetic  field  in  the  form  of  the  superpositioning  of 
oscillations  of  electrical  and  magnetic  types: 
f  .  The  oscillations  of  the  electrical  type 

TM 

are  excited  with  the  currents  ^  ,  and  the  oscillations  of  the 

magnetic  type,  with  the  currents  ,  where  r-r-r  ■ 

Let  us  write  out  the  expressions  for  the  complex  amplitudes 
of  the  component  vectors  of  current  and  field  satisfying  the  in¬ 
homogeneous  Maxwell's  equations: 

xd  H  •  £w„e,£  *J"\  xotE  -  - 1 

divE m  1  ^ 

Using  the  general  expressions  given  in  [13,  2319  we  get  the 
following  relationships  for  one  of  the  spatial  harmonics. 
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Oscillations  of  the  Electrical  Type  (H  =0) 


where 


«-  *8* 

r\  /.i  a/AtO 

c*  TmX  ••  *  ^T5T  * 

£  *  ife/J**  tW«  ' 


The  frequency  of  the  oscillations  is  determined  from  the 


condition 


%(KQ)  •O , 


where  j  n*l,2,3,.... 


In  Eqs.  (2)  and  (4),W,/  is  a  real  function,  finite,  and 
with  its  first  and  second  derivatives  on  the  interval  from 
r=0  to  r=a;  outside  this  interval  it  is  assumed  equal  to  zero; 

P™(cos  8)  are  adjoint  Legendre  functions, 

are  radial  spherical  functions 

r,  e,  <t>  are  spherical  coordinates 
The  coefficients  anm  in  Eqs.  (2)  and  (4)  are  not  defined  here. 
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The  validity  of  the  solutions  presented  is  verified  by  the 
direct  substitution  of  Eqs.  (2),  (3),  and  (4)  in  the  Maxwell's 
equation  (1).  The  vectors  of  the  fields  E*  and  P*  on  the  sur¬ 
face  of  the  sphere  r=a  are  equal  to  zero  and  according  to  the 
theorem  of  equivalency  [13]  the  electromagnetic  field  of  cur¬ 
rents  (2)  in  the  external  region  of  the  system  (r>a)  is  iden¬ 
tically  equal  to  zero. 


Thus,  the  system  of  currents  (2)  is  nonradiating.  The 
power  of  the  electromagnetic  oscillation  generators  is  reactive 
and  can  be  determined  by  the  method  of  induced  electromotive 
forces  with  the  following  formula  [13]: 

P*am^rSJf  £ rf*  (6) 

By  substituting  Eqs.  (2)  and  (3)  in  Eq.  (6)  and  bearing  in 
mind  that 


when  m'+m 
when  m'=m 


J  n(n~<)  P~(c*  •)  P~(cm  • 

f  0  mpm  n’*n 

"  I  when  n'=n, 

we  get  the  following  expression  for  the  generator  power: 


(7) 


We  note  that  the  distribution  of  the  volume  density  of 
electrical  charges  derived  from  the  continuity  equation 

corresponds  to  the  distribution  of  the  volume 
density  of  electric  currents  (2): 

’  -$r  &(*!$*  *' P’(a»  ») 


(8) 


(9) 


Oscillations  of  the  Magnetic  Type  (E  =0) 

J.  -0. 


(10) 

(ii) 

The  frequency  of  oscillations  of  the  magnetic  type  is 
determined  from  the  same  condition  (5).  The  function  fit) 
is  finite  and  continuous  together  with  its  derivatives  on 
the  interval  from  r=0  to  r=a  and  is  equal  to  zero  outside 
this  interval. 

The  validity  of  solutions  (9),  (10),  and  (11)  is  verified 
by  their  direct  substitution  in  Maxwell's  equation  (1).  The 
electromagnetic  field  of  oscillations  of  the  magnetic  type  on 
the  surface  of  a  sphere  r=a  and  outside  this  sphere  is  equal 
to  zero,  that  is,  the  system  of  currents  (9)  is  also  nonradia¬ 
ting. 


where 


H  dp) um.J— 

*  i  ■»»  l  VS St  »  “*  tune  dfdt  • 

f ** —  ML  ft...  dU~ 

•  y/i#  d?  ,  so  * 


SO 


UL  ••Kmt fw  “/WJ t 


Substitution  of  (9)  and  (10)  in  (6)  leads  to  the  following 
expression  for  the  power  of  generators  of  oscillations  of  the 
magnetic  type 


p"« 


(12) 


We  note  that  div j  "»o  ,  and  therefore  the  volume  density 
of  electric  charges  for  oscillations  of  the  magnetic  type  is 
always  equal  to  zero. 
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2.  Superimposing  the  electric  and  magnetic  oscillations 
makes  it  possible  to  reduce  the  system  to  the  resonance  state, 
that  is,  a  state  in  which  the  total  energy  of  the  electromagnetic 
field  of  the  system  with  time  remains  constant  and  changes  from 
electrical  to  magnetic  energy,  and  vice  versa.  This  will  occur 
under  the  condition  that  the  total  reactive  power  of  the  gene¬ 
rators  of  electromagnetic  oscillations  is  equal  to  zero: 

P  -  p**pm.o  (13) 

From  substituting  (7)  and  (12)  in  (13),  we  get  the  expression 


If  we  turn  to  Poynting's  theorem  for  this  system  (as  applied 
to  instantaneous  values) 

'Edtr  •  Trl(*1r~  +  dir,  ( i  5 ) 

then  we  see  that  the  left-hand  side  of  (15),  under  condition 
(14),  proves  to  be  equal  to  zero.  Thus,  the  total  instantaneous 
power  of  generators  of  electromagnetic  oscillations  of  the  sys¬ 
tem  is  also  equal  to  zero  and  the  total  energy  of  the  electro¬ 
magnetic  field  of  the  system  remains  constant  with  time. 

3.  The  functions  /(*)  and  f  (*)  are  arbitrary  in  the  sense 
indicated  above.  We  can,  in  particular,  assume  t  . 

Then  condition  (14)  can  be  reduced  to  the  following: 


If  the  following  boundary  conditions  are  imposed  on  these 
functions : 


•  0 


IX 


(16) 
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the  electrical  charges  and  currents  at  the  boundary  of  the 
system  (r=a)  prove  to  be  equal  to  zero.  So  they  will  be  con¬ 
tinuous  everywhere,  including  the  boundary  of  the  system. 
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ITERATIVE  METHODS  OF  SYNTHESIZING  REFLECTION  ARRAYS 
A.  F.  Chaplin,  V.  A.  Mashkov,  and  V.  F.  Makkaveyeva 

This  paper  examines  certain  iterative  methods  of  synthe¬ 
sizing  methods  of  synthesizing  arrays  of  radiators  leaded 
with  reactances.  Because  there  is  no  widely  accepted  name 
for  this  class  of  arrays,  we  will  call  them  "reflection" 
arrays,  as  is  done  in  a  number  of  publications  .  The  formu¬ 
lation  of  the  problem  of  synthesizing  these  arrays  was  given 
in  the  study  [25].  In  examining  reflection  arrays,  the  system 
of  Kirchhoff  equations  is  conveniently  written  out  in  the  form 

{CZJ  *  CE3  <X>}  <t>  ■  <r>  ,  (1) 

where  [Z]  is  the  matrix  of  the  assigned  Intrinsic  and  muutal 
impedances  of  the  radiators 
[E]  is  the  unit  matrix 

<X>  is  the  vector-column  of  the  load  reactances 
<i>  is  the  vector-column  of  currents  in  the  radiators 
<e>  is  the  vector-column  of  the  voltages  induced  by 
the  active  radiator. 

The  solution  of  this  problem  of  synthesis  can  be  reduced 
to  finding  the  reactances  <X>  based  on  the  given  radiation 
pattern  of  an  array;  here  the  currents  <i>  must  belong  to  the 
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set  J  providing  the  reactance  of  the  impedances  <X>.  The  appro¬ 
ximation  to  the  given  radiation  pattern  can  be  root-mean-square 
of  uniform.  The  problem  of  the  root-mean-square  approximation 
for  a  linear  array  consists  of  minimizing  the  function 

ru„...  o  •/  lr(»)  -  9(9)  Z  im t**'**lid9h  ( 2 ) 

and  the  uniform  approximation  problem,  to  minimizing  the 
function 

n fi«,  i„...  O  •  moxlF(9)-C9(9)LLmtfm4“*U,9\  .  . 

9  *•'  ( 3 ) 

The  final  formulation  of  the  problem  of  synthesizing  a  reflec¬ 
tion  array  with  arbitrary  fixed  geometry  can  be  represented 
as  follows:  minimize  T  (1^,  i^ >  ...»  i-n)  with  respect  to  <i> 
under  the  condition  that  <i>ej  >  that  is,  find  the  <3f>  such 
that 

r(Z.  V •  ••  O  ■  yis  T(i„  ...  L J, }  (H) 

<Z>€  J  I 

Ii.  mathematics  there  is  not  yet  an  algorithm  capable  of 
solving  problem  (4)  in  the  most  general  form,  that  is,  when 
r(ilt  .  ..,  in)  is  an  arbitrary  function,  and  J  is  a  set  with 
arbitrary  configuration.  In  the  case  of  reflection  arrays, 
the  currents  from  the  J  must  satisfy  the  following  system  of 
constraints  [25]: 


m  V  CJ  m°,  *•<.  N 


(5) 


This  article  presents  the  results  obtained  using  different 
algorithms  developed  by  the  authors  as  applied  to  the  problem 
of  reflection  arrays. 
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Let  us  assume  we  know  certain  initial  values  of  the  load  im- 
pendances  of  the  dipoles  <x>^°\  Selection  of  the  initial  value 
of  the  loads  can  be  done  with  different  methods: 


a)  First  Method  of  Selecting  Initial  Loads 

Let  us  solve  problem  (4)  without  allowing  for  the  constraints 
on  the  currents.  For  this  purpose  it  is  best  to  set  up  an  over¬ 
determined  system  of  equations  that  can  be  solved  for  function  (2) 
by  the  method  of  least  squares,  and  for  function  (3),  by  the 
Stieffel  exchange  method  [24].  After  finding  the  currents  <i>, 
which  generally  do  not  belong  to  the  set  J ,  let  us  calculate  the 
load  impedances  based  on  a  linewise  solution  of  system  (1)  and, 
by  eliminating  the  real  part  in  this  solution,  we  regard  them  as 
the  initial  loads. 


b)  Second  Method  of  Selecting  Initial  Loads 


Conditionally  we  assume  that  the  mutual  influence  between  the 
dipoles  in  the  array  is  weak,  that  is,  the  matrix  [Z]  is  a  diagonal 
matrix.  In  this  case,  the  currents  from  (1)  are  defined  as 


■ 


(6) 


The  right-hand  side  of  Eq.  (6)  is  a  circle  on  a  complex  plane 
[25].  If  we  denote  the  current,  which  is  the  solution  of  problem 
(4)  without  the  presence  of  constraints  (5),  by  <i*>,  as  the  ini¬ 
tial  value  of  the  load  impedances  we  select  the  value  that  mini¬ 
mizes  the  expression 


K  • 


(7) 


The  solution  of  Eq.  (7)  has  a  simple  geometrical  significance 
and  analytically  can  be  reduced  to  solving  a  quadratic  equation. 
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The  simplest  algorithm  for  solving  the  problem  of  synthe¬ 
sizing  an  array  in  the  root-mean-square  approximation  is  given 
in  [26].  The  problem  can  be  reduced  to  minimizing  the  function 
A: 

-A,  (8) 

given  the  system  of  constraints  on  current  (5).  The  algorithm 
above  presupposes  that  the  matrix  [Z]  has  an  emphasized  diagonal 
and  the  current-column  <i>^0^  is  known,  which  we  adopt  as  the 
zero  approximation.  The  assumption  that  there  is  weak  mutual 
influence  between  the  dipoles  enables  us  to  minimize  (8)  for 
the  k-th  dipole,  when  there  is  the  k-th  constraint  on  system  (5). 
The  Lagrange  method  as  applied  to  this  problem  leads  to  a  system 
of  two  equations,  one  being  quadratic,  and  the  other,  linear. 

The  solution  obtained  in  this  case  does  not  exactly  satisfy 
conditions  (8)  and  (5).  For  a  revised  solution,  new  computa¬ 
tional  cycles  are  carried  out.  This  algorithm  makes  it  possible 
to  compute  the  arrays  of  a  small  number  of  dipoles  quite  widely 
dispersed  in  space  (d/A>0.1).  Its  advantage  is  the  simplicity  of 
computations.  The  desire  to  allow  for  the  influence  of  only  two 
adjoining  dipoles  leads  to  much  more  complexity  of  computations. 
Elaborating  this  concept,  we  now  must  minimize  (8)  for  the  k-th 
dipole  given  the  (k-l)-th,  k-th,  and  (k+l)-th  constraints  on  the 
current  from  system  (5).  The  Lagrange  method  leads  to  a  system 
of  nine  equations,  three  of  which  are  quadratic.  An  incomplete 
allowance  for  the  mutual  influence  leads  to  a  situation  in  which 
the  solution  obtained  will  be  inexact;  and  just  as  in  the  first 
case,  repeated  computations  will  be  necessary.  When  the  mutual 
influence  between  all  dipoles  is  taken  into  account,  a  system  of 
3N  equations  must  be  solved,  of  which  N  are  quadratic.  In  con¬ 
trast  to  the  first  two  cases  examined,  with  an  increase  in  the 
number  of  elements,  the  matrix  order  also  grows  larger.  Allow¬ 
ing  for  mutual  influence  (partially  or  completely)  leads  to  the 
need  to  solve  a  system  of  nonlinear  equations.  Existing  mathe¬ 
matical  methods  enable  this  problem  to  be  solved  [27]. 
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We  present  our  further  exposition  on  the  example  of  the  uni¬ 
form  synthesis  of  a  reflection  array.  For  the  case  of  the  root- 
mean-square  synthesis,  the  algorithm  functions  quite  analogously; 
there  will  be  a  difference  only  in  the  normalizing  of  the  radia¬ 
tion  pattern.  The  problem  of  uniform  synthesis  is  formulated  in 
the  form  of  problem  (4),  where  r(i1,...,in)  stands  for  function 
(3)>  under  the  condition  that  the  set  D  is  defined  by  (5)-  Into 
function  (3)  is  introduced  the  complex  normalizing  coefficient  C, 
which  must  provide  us  with  the  convergence  of  the  sequence  <i>^0^ 
<i>^^,  ...,  <in>  t0  the  desired  distribution  <£>€J  .  Physically 
it  signifies  the  possibility  of  adjusting  the  radiation  of  the 
array  to  the  level  that  is  optimal  in  the  sense  of  the  uniform 
proximity  of  two  curves:  the  normalized  given  radiation  pattern 
F  (  0 )  : 

max  F{*)*  1 ;  •A»C"5,*^3  (9) 

and  the  radiation  diagram  of  the  array  f(0).  The  magnitude  of 
the  coefficient  C  Is  selected  by  solving  the  following  auxiliary 
problem  for  a  known  (calculated)  current 

4  •  maxima) I  * 

•  mi**  max 
C  ^ 

The  solution  of  problem  (10)  Is  obvious: 

?m  rtf'?*  + 

♦  Umax 

The  soundness  of  the  introduction  of  this  coefficient  follows  from 
the  fact  that  the  left-  and  right-hand  members  of  Eq.  (1)  can  be 
multiplied  by  an  arbitrary  number.  The  essentials  of  the  proposed 
algorithm  will  consist  in  constructing  a  sequence  of  current  dis¬ 
tributions,  each  term  of  which  has  a  smaller  quantity 


(10) 


(11) 
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i".  *».m- i . 

0 

( k_l ) 

compared  with  A 

To  calculate  the  quantity  A,  we  partition  the  interval  of 
angles  on  which  F(9)  and  f(9)  is  defined  with  a  closely  spaced 
net  and  we  calculate  A  at  the  nodes  of  this  net.  Let  us  adopt 
the  calculated  initial  values  of  load  impedances  */**/*.. X.J* 
and  insert  them  into  Eq.  (1).  From  the  system  of  equations  (1) 
we  calculate  the  current 

(12) 

From  Eq.  (11)  we  calculate  the  normalized  coefficient  cor¬ 
responding  to  these  currents,  and  then  the  quantity  In 

the  following  text  we  leave  unchanged  all  the  initial  values  of 
the  load  impedances  of  the  dipoles,  except  for  the  load  impedance 
of  the  first  dipole.  By  varying  the  quantity  jX,  ,  each  time  we 
calculate  the  values  of  <i>,  ft,  and  A,  attempting  to  find  a  » 
such  that 

(13) 

To  calculate  the  .currents  corresponding  to  each  new  value  of  the 
quantity  jX,  generally  speaking  each  time  we  must  invert  the  matrix 

[LZJ+ttXty  . 

If  the  number  of  attempts  for  the  search  will  be 
large  and  the  matrix  is  large  in  size,  this  procedure  requires 
more  machine  time.  This  difficulty  is  cleared  away  with  two 
procedures : 

1)  The  matrix  {CZW£XX>}  corresponding  to  each  new  value  of 
can  be  obtained  not  by  inversion,  but  by  correcting  the 
Inverted  matrix  corresponding  to  the  previous  value  of  yX/""*. 

The  correction  is  carried  out  using  the  following  formulas  [23]. 
Let 
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and 


.  {cz]*C£j<x>*T-r«3 
{rzj  -  [£]  <x4Mx-_  xr>4 f.  cmj  “* 


(14) 


t  hen  mr%~  £/TVi  /r,/» 

where  £  •  7$^  ,d‘j(X?-X?) 


2)  The  strategy  of  searching  for  must  be  such  that  find¬ 
ing,/^  occurs  in  the  smallest  number  of  attempts.  This  optimal 
strategy  can  be  developed  only  by  relying  on  the  form  of  the 
function  *(/*,)  .  For  an  arbitrary  number  of  dipoles,  the  form 
of  this  function  cannot  be  analytically  investigated,  but  for 
two  dipoles  it  is  shown  in  Fig.  1.  For  the  form  of  function  shown 
the  optimal  strategy  was  constructed  using  a  Fibonacci  series. 

The  function  was  computed  at  the  following  points: 


fX",  j  (X?*  {««)#  •  •  • 

c,  a;4. 


ZBo) ; 


(15) 


In  Hq.  (15)  is  the  k-th  member  of  the  Fibonacci  series  1,  1,  2, 
3,  5,  8,  . . . ;  a  is  the  initial  interval,  equal  to  +0.05;  ±0.5;  and 
+5,  which  is  derived  in  order  to  occupy  a  section  to  the  left  of 
.  The  condition 


(16) 

is  a  criterion  for  converting  from  one  term  of  the  sequence^fc? 
another.  The  search  for  ends  when  (16)  is  not  satisfied  for 
any  values  of  the  parameter  a=+0.'05,  +0.5,  and  +5. 


To  continue  the  process,  we  now  record  the  value  of  found 
and  carry  out  an  analogous  procedure  for  ,  where  all  jX™ 

remain  unchanged,  as  before.  All  the  load  Impedances  are  selected 
in  this  way. 

We  have  described  one  cycle  of  the  algorithm.  The  entire 
algorithm  consists  of  these  repeating  cycles.  A  quantity 
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analogous  to  ( 1 6 )  is  the  criterion  for  termination.  The  conver¬ 
gence  of  this  process  to  <1>  can  be  easily  demonstrated  if  v*. 
has  the  form  shown  in  Pig.  1. 


As  an  example,  let  us  examine  the  synthesis  of  an  equi¬ 
distant  reflection  array  if  the  given  radiation  pattern  has  the 


following  form: 

I  i,  -?***¥ 


(17) 


Considering  the  plane  0  as  the  plane  of  isotropicity  of  the 
lattice  *4-1,  we  get  the  following  results  by  applying  the  above- 
described  algorithm  to  problem  (3)-(^)*  For  N=2$  and  d=0.25^, 
when  the  number  of  points  of  partitioning  of  the  interval 
[0,tt/2]  is  equal  to  m=19,  the  desired  load  impedances  <X>  in 
ohms  as  shown  in  Table  1  are  obtained.  The  assigned  and  the 
obtained  radiation  patterns  are  presented  in  Pig.  2.  Prom  Fig.  2 
we  see  that  the  approximation  proves  to  be  quite  good  if  we  omit 
consideration  of  the  large  side  lobe. 


In  conclusion  we  must  note  that  all  the  algorithms  proposed 
in  this  article  can  be  classed  with  local  methods  of  analysis  of 


multiextremal  problems.  To  be  fully  convinced  that  the  minima 
found  are  global,  we  must  apply  more  complex  nonlocal  methods 
of  analysis,  for  example,  the  I.  M.  Gel’fand  method  of  troughs. 
Naturally,  converting  to  nonlocal  methods  is  necessary  when  the 
magnitude  of  the  local  minimum  does  not  suit  the  investigator. 

Table  1 

Xii  |Xtf  |x±8  ?X*  fey  \x&  jx*7  jy*a  jy±9'.y^o 

-;22  l->8  l->a  lyao  I  y  7, 31-/1037 18  I-/  5  lyisiyia 
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ENERGY  METHOD  OF  SYNTHESIZING  MICROWAVE  ANTENNAS 
A.  F.  Chaplin 

At  the  present  time,  the  problem  of  finding  the  current 
distribution  or  the  distribution  of  the  tangential  field  com¬ 
ponent  in  the  antenna  aperture  for  a  given  radiation  pattern 
has  been  formulated  and  solved  in  the  literature  on  antenna 
synthesis.  It  must  be  stated  that  with  this  formulation  of  the 
problem,  the  antenna  must  be  represented  as  the  distribution  of 
extrinsic  current  (electric  or  magnetic)  sustained  by  current 
or  voltage  generators.  At  the  same  time  we  know  that  micro- 
wave  antennas  are  characterized  generally  by  a  free  excitation 
mode  [29],  in  which  the  antenna  is  fed  a  flow  of  incident  power 
arriving  from  a  constant  power  output  generator.  In  addition, 
with  the  usual  formulation  of  the  antenna  synthesis  problem, 
estimates  of  antenna  energy  characteristics  grow  in  complexity. 
This  is  evidenced  by  the  fact  that  in  most  studies  on  synthesis 
the  antenna  power  is  calculated  as  the  integral  of  the  square  of 
the  radiation  pattern  module  with  respect  to  field  intensity. 

We  can  easily  show  that  this  estimate  is  invalid  above  all  in  the 
calculation  of  the  total  and  reactive  antenna  power.  Examined 
below  is  a  formulation  of  the  antenna  synthesis  problem  in  which 
there  is  an  association  between  the  distribution  of  normalized 
voltage  in  the  far  zone  and  in  the  antenna  aperture. 
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As  a  working  model,  we  adopt  a  linear-array  antenna  with 
length  21,  arranged  along  the  y-axis  and  irradiating  a  field  in 
the  form  of  TM-waves  with  respect  to  the  z-axis  perpendicular 
to  the  antenna.  Let  us  examine  the  two-dimensional  problem  when 
the  field  does  not  depend  on  the  x  coordinate.  The  free-space 
field  is  represented  in  the  form  of  a  superpositioning  of  plane 
waves;  to  do  this,  we  use  expressions  from  the  monograph  [13]: 


(1) 


Here  it  is  assumed  that  the  observation  point  is  located  above 
the  extrinsic  sources,  whose  spectral  density  is  represented 
by  the  function  F(*)- 

Each  plane  wave  produces  a  power  flow  perpendicular  to  the 
antenna  and  equal  to 


Let  us  represent  the  quantity  II  in  the  form  of  the  product  of  two 
symmetrical  cofactors 

fi,  "  f  *(»)  fy*)  , 

where 

ore  ■'far  ;  ow-flsr  (2) 

This  representation  is  similar  to  the  representation  widely  used 

in  microwave  circuit  theory  for  normalized  voltage  waves  in  wave- 

1/2 

guides.  The  quantity  0(«)  has  the  dimension  of  (watts) 

As  we  know,  the  current  distribution  along  the  linear-array 
antenna  is  associated  with  the  field  intensity  in  the  far  zone  by 
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means  of  a  Fourier  transform.  In  the  same  way,  we  can  associate 
the  function  $>(y)  and  4(a),  which  are  the  power  distributions  to 
the  exponent  1/2  along  the  linear-array  antenna  and  in  the  far  zone, 
respectively . 


4(a)  •  f<P(y)  tim*  dy, 
it  — 

d* 

If  the  antenna  is  located  in  loss-free  space,  the  values  of  the 
function  *(*)  along  the  segment  [-k,k]  define  the  active  radiation 
power,  and  the  values  of  this  function  at  the  semiinfinite  segments 
[-°°,  -k]  and  [k,°°]  define  the  reactive  energy  stored  in  the  near 
field  of  the  antenna. 

Thus,  the  problem  of  synthesizing  microwave  antennas  can  be 
posed  and  solved  as  follows:  we  assign  the  function  9(a)  and  with 
the  solution  of  the  integral  Fredholm  equation  of  the  first  kind 
(3a)  we  find  the  function  <t(y)  from  which  we  can  calculate  the 
distribution  of  the  flow  of  complex  power  along  the  antenna.  If 
the  function  9(a)-  is  defined  along  the  entire  «  axis,  to  solve  the 
synthesis  problem  we  can  use  the  inversion  formula  (3b). 

To  solve  the  problem  posed  we  can  develop  methods  similar 
to  methods  of  solving  the  conventional  antenna  synthesis  problem: 
the  Fourier  integral  method,  the  method  of  partial  radiation  pat¬ 
terns,  and  the  A.  N.  Tikhonov  regularization  method.  It  must  be 
noted  that  now  all  energy  estimates  will  have  an  explicit  physical 
meaning. 

In  addition  to  calculating  the  flow  of  complex  power  along 
the  antenna,  we  can  find  the  distribution  of  the  surface  impedance 
or  the  input  impedance  at  each  antenna  point.  The  surface  impedance 


(3a) 

(3b) 
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Z(y)  can  be  found  with  the  formula 

f,NI-  < 

Let  us  dwell  cn  two  properties  of  the  functions  W*)  and 
(y)  we  have  introduced 

1.  Function  *(■)  and  My)  exhibit  the  property  of  local 
separation  of  power  into  active  and  reactive  at  each  antenna 
point.  Let  us  take  two  plane  waves  with  the  wave  numbers 

*»  and  xt  :  +,(*}- At (k-k,),  .  Let  and /*,/>*. 

Let  us  find  the  conditional  power  flow  at  some  point  p(y,0)  of 
“he  antenna. 


tAliaicosCax 9  A-argg 
•Ik1-** 


The  power  flow  of  two  slow  plane  waves  represented  with  the  func¬ 
tions  4>^(y)  and  remains  reactive  at  each  point  on  the  plane 
z=0 . 


Let  iKh*  and  faej<*  •  Then  we  have  the  expressions 

♦  IM1SL  cosCaiaA  -atgB 

The  power  flow  of  two  fast  plane  waves  represented  by  the  functions 
<{> 2 ( y )  and  $2(y)  remains  active  at  each  point  on  the  plane  z=0. 

3ut  if  we  examine  the  power  flow  of  one  fast  and  one  slow  wave 
and  ,  it  proves  to  be  complex  at  some  point  p(y,0): 
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2.  The  transition  to  the  normalized  voltage  waves  4>(y)  signi¬ 
fies  the  symmetrizing  of  the  kernel  of  the  integral  equation,  which 
associates  E  (y)  (or  the  surface  magnetic  flux  3* (V))  with  the  dis¬ 
tribution  of  the  surface  impedance  Z(y)  and  can  be  derived  from 
Eq.  (4): 

mm 

£s  (v)  Z(y j  dv'  ( 5 ) 

Actually,  by  designating  mCPN  >  we  bring  Eq.  (7)  to 

the  form  — 

&(*)  -t4-  J {Z(?)'iZMH*>(i<iy-vl)<P(v’)d'/' , 

where  the  kernel  is  already  symmetrical. 
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DIFFRACTION  OF  PLANE  WAVES  BY  AN  ARBITRARY  SMOOTH  CONVEX  BODY 
AND  THE  KELLER  METHOD  OF  DIFFRACTION  WAVES 

N.  S.  Orlova  and  Yu.  I.  Orlov 


A  shortwave  asymptotic  representation  for  an  electromagnetic 
field  was  derived  in  [30],  which  is  valid  in  the  boundary  layer 
near  an  ideally  conducting  surface  S  in  the  region  of  deep  shadow 
and  semishadow.  Using  the  expressions  for  currents  at  the  sur¬ 
face  of  a  body  derived  in  [30],  in  this  paper  the  diffraction 
field  in  the  wave  zone  in  the  shadow  region  is  determined.  The 
solution  obtained  generalizes  the  Keller-  formula  [31]  for  the 
region  of  shadow  and  is  utilized  for  refining  the  concept  of 
diffraction  rays  (the  geometrical  theory  of  diffraction)  [31]. 
Analogous  problems  in  the  two-dimensional  case  were  investigated, 
in  particular,  in  the  studies  [33-37]. 

1.  Scalar  Problem 


Suppose  that  at  a  smooth  convex  surface  a  plane  wave  is 
incident  at  .  Let  us  examine  the  first  and  second  boun¬ 

dary  value  problems  for  the  Helmholtz  equation 


U  •  L /'♦  U\  §ff'  L*Um.  o(r«)t 
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where  a=l,  for  the  Dirichlet  problem,  and  St ■ ^  ,  for  the  Neumann 
problem  (n  is  the  normal  to  the  surface  S). 

Near  the  surface  S  the  field  U(P')  in  the  shadow  region  can 
be  defined  with  the  formulas  given  in  [30].  In  this  case  we  use 
the  coordinates  of  the  boundary  layer  v,u,n  (Pig.  1),  where  u= 

=  QO.*Q.P.,  the  shortest  distance  between  point  Q  on  the  incident 
wavefront  and  the  point  Pq,  and  v  •  Q~Q9  ,  the  distance  along  the 
boundary  of  the  geometrical  shadow  (Q-^  is  the  initial  reference 
point),  and  n=P'Pg,  the  distance  along  the  normal  n.  (  w  »•*)  to 
the  surface  S,  whose  equation  is  r*r%(u  ,tr)  ,  The  system  of  coor¬ 
dinates  u,  v  on  S  is  semigeodetic ;  in  this  case,  the  first  prin¬ 
cipal  quadratic  form  of  the  surface  is  of  the  form 
d3‘»du?.hl(u,*)<**4  >  where  K,  C  )'  -  I. 

To  define  the  field  in  the  shadow  region  far  removed  from 
the  surface  let  us  use  Green's  formula  [38],  with  can  be  written 
approximately  as  follows,  with  reference  to  the  results  in  [39]: 


where  R»lr-r?(u,*)l  ,  and  the  integration  is  performed  over  the 

shadow  part  S,_  of  surface  S.  Here  and  everywhere  in  the  following 
1/ 

text  the  first  (second)  row  in  the  double  formulas  applies  to  the 
Dirichlet  (Neumann)  problem. 

The  expressions  for  the  quantities  and  Ug ,  which  appear 
in  Eqs.  (2)  can  can  be  called  the  "currents"  I,  follow  from  the 
boundary  layer  formulas  given  in  [30]  when  m=0: 


where  f#(Z)  and  g*(Z)  are  functions  complexly  adjoint  to  the 
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Foch  functions  [40] 


P  is  the  radius  of  curvature  of  a  geodetic  (ray),  and  v=const; 
in  Eq.  (3)  and  in  the  following  treatment 


V*.®  J“J>W  and  %.u(Q.) 

^"^^(Fig-  D*  The  region  of  applicability  for  Eqs .  (3)  for 
currents  is  the  region  of  shadow  and  semishadow,  all  the  way  to 
merger  with  geometrical  optics  formulas.  By  inserting  (3)  into 
(2),  when  kR>>l  we  get 


(4) 


where  I'-RR  ,  and  A  is  the  region  of  variation  of  u  and  v  corres¬ 
ponding  to  the  shadow  part  of  the  surface  ST. 

The  double  integral  in  (4)  we  calculate  with  the  stationary 
phase  method.  Here  the  classical  stationary  phase  method  [39] 
turns  out  to  be  inapplicable.  So  first  the  interior  integral  in 
v  is  computed  asymptotically,  and  then  the  refined  stationary  phase 
method  is  applied  to  the  resulting  integral  in  u  [37].  As  a  re¬ 
sult  we  get  the  following  expressions  for  the  field  in  the  shadow 
region  at  large  enough  distances  from  the  surface  S  (kRg>>l). 


A  £  AJJ 

where  f  (ZQ)  and  g  (ZQ)  are  functions  complexly  adjoint  to  the 
Foch  functions  [37»40,4i] 

R0  Is  the  length  of  the  tangent  PPQ, 
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u,,  UJ  are  the  coordinates  of  the  separation  point  of  the 
diffraction  ray  . 

&D 

Eqs.  (5)  are  a  generalization  of  Keller's  geometrical  dif- 

A 

fraction  formulas,  since  they  are  valid  also  in  the  semishadow 
region. 

Vector  Problem 

Let  us  examine  the  problem  of  the  incidence  of  a  plane  elec¬ 
tromagnetic  wave  by  an  arbitrary,  smooth 

convex,  ideally  conducting  body.  The  field  in  a  boundary  layer 
with  thickness  xn&H  »  I)  near  the  surface  of  the  body 

in  the  shadow  region  is  determined  by  the  equations  in  [30].  To 
determine  the  field  in  the  shadow  region  far  from  the  body  sur¬ 
face  let  us  use  the  equivalency  theorem  [13],  which  with  refer¬ 
ence  to  [39]  can  be  approximately  represented,  analogously  to 
(2) ,  in  the  form: 

T=r  ( 6  > 

'  (7) 


where  I-fsjwJ/,  is  the  surface  density  of  the  electric  current  in 
the  shadow  part  ST  of  the  metal  surface  S.  In  accordance  with 
[30]  we  have: 


(8) 


Here  jr.  |g.  and  are  the  unit  vectors  of  the  semigeodetic 

system  of  coordinates  u,v,  where  Ej*-  *  are  the 

components  of  the  incident  field  at  the  boundary  points  QQ 
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(Pig.  1),  Hipjj  and  are  the  "currents"  defined  by  Eqs.  (5)  in 

the  Dirichlet  and  Neumann  problems,  respectively. 


Prom  Eqs.  (6)  and  (7),  with  reference  to  Eq.  (8),  as 
kR>>l,  we  get 


Ca.  -j-  dS, 


(9) 


where  O,- ST".  •  By  amoving  the 

slowly  varying  cofactors  from  under  the  integral  sign  at  the 
stationary  point  u0,vQ,  defined  by  the  condition  ,  from 

Eq.  (9)  with  reference  to  (2)  it  is  not  hard  to  find  the  follow¬ 
ing  expressions  for  the  field  in  the  shadow  region  at  large 
enough  distances  from  the  surface  S  (kRQ>>l) 

I(r)'£^U,(nn,(u„%)  ♦  E'^CiPWu..  *.), 

H (F)  •ffiCell(u**,)E(r)]  •  ( 10 ) 

•  HHlkt/4(r)n,(u^  t*)  *  (F)  *<), 

where  are  the  normal  components  of  the  inci¬ 

dent  field  at  the  boundary  points  QQ  (Pig.  1),  U^( r)  and  U2(r)  of 
solution  (5)  of  the  Dirichlet  and  Neumann  problems,  n.(u* VJ  and 
Tjpj^vQ  are  the  unit  vectors  of  the  principal  normal  and  of  the 
binormal  of  the  geodetic  \-cantt  at  the  separation  point  PQ  of 
the  diffraction  ray  PPQ. 

The  region  of  applicability  of  Eqs.  (10)  is  the  same  as  in 
the  scalar  case. 
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DIFFRACTION  OF  A  PLANE  ELECTROMAGNETIC  WAVE  BY  AN  IDEALLY 
CONDUCTING  BODY  OF  REVOLUTION 

N.  S.  Orlova 


1.  Let  the  plane  electromagnetic  wave 

(1) 

be  incident  on  a  smooth,  convex  surface  of  revolution  S,  whose 
radii  of  curvature  are  much  larger  than  the  wavelength,  and  the 
equation  of  the  generatrix  in  a  cylindrical  system  of  coordinates 
is  of  the  form 


or  '*•*'<’*)  (2) 

According  to  [30],  let  us  Introduce  near  S  the  semigeodetic 
system  of  coordinates  v,u,n  normally  associated  with  it  and 
formed  by  the  single-parametric  family  of  geodetlcs  v=const 
and  the  curves  u=const  orthogonal  to  It;  n  is  the  length  of  the 
segment  of  the  normal  to  the  surface.  As  we  know  [42],  the  tra¬ 
jectories  of  the  geodetic  lines  at  a  convex  body  of  revolution 
fluctuate  between  two  parallels  r=h  and  are  governed  by  the 
following  equations 


rtinji  "h, 

df  .  h  iL 

3umP,  Su  r 


(3) 

(4) 
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where  du  is  an  element  of  the  length  of  the  geodetic,  8  is  the 
angle  between  the  geodetic  and  the  meridian,  and  h  is  the  radius 
of  the  smallest  parallel  that  is  tangent  to  the  given  geodetic. 

The  geodetic  lines  v=const  are  envelope  rays  at  the  surface 
issuing  from  the  boundary  points  of  the  geodetic  shadow  on  S 
and  serves  as  an  extension  of  the  incident  rays 

(5) 

(t^  is  the  unit  vector  of  the  envelope  ray). 

The  equation  of  the  boundary  of  the  geometric  shadow  at  the 
surface  is  found  by  determining  the  common  points  on  the  surface 
S  and  the  family  of  planes  tangent  to  it  and  parallel  to  the 
vector  of  the  wave  normal  of  the  incident  wave  ic 

(6) 

where  V  #!«.«* 

By  solving  jointly  Eqs.  (3)  and  (5),  let  us  determine  the 
parameter  of  the  geodetic  #*•*(■*  J?  ) 

IN  ^  ^  ^ 

Thus,  for  the  geodetlcs  v=const,  Eqs.  (4)  can  be  written  in 
the  form 


(8) 

(9) 

where  9^  and  are  determined  by  Eq.  (6),  and  is  the  phase 
of  the  incident  field  at  the  exit  point  of  the  geodetic  at  the 
boundary  of  the  shadow 


11  jsm 


Eqs.  (8)  and  (9)  are  valid  when  where  the  angle  of  rota¬ 
tion  corresponds  to  the  tangency  point  of  the  geodetic  with 

the  parallel  r=h.  When  ,  in  Eqs.  (8)  and  (9)  we  must  use 

the  partitioning  of  the  integrals  in  terms  of  the  rotation  point, 
with  reference  to  the  symmetry  of  the  geodetic  with  respect  to 
. 


As  the  coordinate  v  it  is  convenient  to  select 

(ID 

It  is  not  difficult  to  show  that  the  curvilinear  coordinates 
v  and  u  selected  are  orthogonal.  The  system  of  coordinates 
v,u,n,  generally  speaking,  is  orthogonal  only  near  the  surface. 
The  Lame  coefficients  of  the  resulting  semigeodetic  coordinate 
system,  for  small  n,  can  be  expressed  as  follows: 


where 


.  • 


(12) 


2.  Following  the  study  [80],  let  us  write  out  expressions 
for  the  diffraction  fields  f  and  S  near  the  surface 

{«: 


Here 


(13) 


R  R 

J>  —  is  the  radius  of  curvature  of 

the  geodetic,  R1  and  R2  are  the  principal  radii  of  curvature  S 


71 


y/fe K 9]  are  functions  complexly  adjoint  to  the  Foch  func¬ 
tions  [43] 


The  quantities  and  ^  can  be  calculated  at  the  exit  point  of 
the  diffraction  ray  from  the  shadow  boundary.  £  and  are 
determined  by  the  expansion  of  the  components  of  the  incident 
field  (1)  with  respect  to  the  coordinates  v,u,n  at  the  shadow 
boundary 


where 


-«-•  -i-fj  ■  % 
4*  itn<e*eoixA^ 


(14) 


a*  • 


The  conversion  from  the  semigeodetic  system  of  the  coordi¬ 
nates  v,u,n  to  the  spherical  coordinates  9,$  at  the  surface  is 
carried  out  based  on  the  formulas 


(15) 


Eqs.  (13)  define  the  field  corresponding  to  a  single  dif¬ 
fraction  ray  arriving  at  the  observation  point.  In  the  simplest 
case,  the  field  in  the  shadow  region  of  the  body  of  revolution 
is  formed  by  means  of  two  diffraction  rays  enveloping  the  surface 
on  opposite  sides.  But  the  ray  pattern  of  the  field  depends  on 
the  specific  form  of  the  surface  and  in  the  general  case  can  be 
quite  complex.  In  particular,  investigation  of  the  pattern  of 
diffraction  rays  at  the  ellipsoid  of  revolution  conducted  as 
per  Eq.  (8)  with  a  numerical  method  shows  that  the  number  of  rays 
arriving  at  the  arbitrary  point  P  of  the  surface  depends  strongly 
on  the  elongation  of  the  ellipsoid  and  on  the  angle  of  incidence 
of  the  plane  wave  (Fig.  1). 


Thus,  calculation  of  the  fields  at  each  specific  surface  must 
be  preceded  by  the  determination  of  the  number  of  diffraction 
rays  arriving  at  different  zones  of  the  shadow  region.  The  fields 
of  all  rays  are  determined  by  Eqs.  (13)  and  differ  from  each  other 
only  by  the  values  of  the  variables  u  and  v  arriving  there. 
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t-*r 


-  wi* 


.  -  shadow  boundary 

Pig.  1 


Fig.  2 


3 


3.  The  solution  (13)  that  we  derived  cannot  be  used  near 
the  caustics  (when  h  ->-0).  In  the  cases  when  the  ray  pattern  at 
the  surface  has  a  caustic,  the  solution  (13)  must  be  analytically 
extended  into  the  caustic  region,  for  example,  using  the  method 
given  in  [44]. 


In  the  case  of  the  axial  incidence  (9°=0)  of  a  plane  wave, 
an  axial  caustic  always  originates  on  the  body  of  revolution  of 
arbitrary  shape;  to  determine  the  field  near  the  axial  caustic 
Eqs.  (13)  can  be  similarly  modified  [45].  In  this  case,  h=0  and 
the  diffraction  rays  are  propagated  along  the  meridians,  so  it  , 
is  sufficient  to  examine  the  ray  pattern  in  the  plane  9>9r»cons4 
(Pig.  2).  At  each  observation  point  two  rays  arrive,  each  pro¬ 
pagating  along  the  shortest  path,  and  the  other  envelopes  the 
"dark  pole"  (z*-zm_„),  thereby  acquiring  an  additional  lead  in 
phase  ir/2. 


As  we  can  see  from  Pig.  2  and  Eq.  (12),  for  both  rays  we 
have  the  expressions 

h*mir.  r~  ‘ (16) 


The  coordinates  of  the  rays  arriving  at  point  P  from  points  F-j-  and 
F-j-j  can  be  written  in  the  form 


(17) 


where  '«■*  is  the  ray  coordinate  of  the  rotation  point  z=-z 

r  mzx 


Assuming  that  the  following  is  true  near  the  axial  caustic 

and  referring  to  Eqs.  (16)  and  (17)  and  the  additional  lead  in 
phase  at  the  caustic,  let  us  represent  the  composite  field  near 
the  axial  caustic  as  follows: 


74 


c 


(13) 


(19) 


Using  the  familiar  expansion  of  Bessel  functions  [46] 

3jx)  •  {fx  LCOS(x-  T  *  fj  *  0&)]j 

from  Eq.  ( 1 8 )  we  get 


w)  •  <20) 

Here  and  are  calculated  at  any  point  on  the  shadow  boun¬ 
dary.  We  note  that  . 

Similarly,  by  referring  to  the  mutual  orientation  of  the 
vectors  v,u,n  for  the  first  and  second  rays,  we  get 

,y)’ 


We  note  that  Eqs .  (20)  and  (21),  in  contrast  to  Eqs.  (13), 
describe  the  total  field,  considering  all  the  rays  focussing 
on  the  axial  caustic. 

The  expressions  we  have  derived  for  the  diffraction  field 
of  the  plane  wave  are  valid  near  the  surface  of  the  body  of 
revolution  with  an  arbitrary  shape;  in  the  case  of  a  sphere 
they  wholly  coincide  with  the  asymptotic  exact  solution  [47], 
Including  the  caustic  region.  The  resulting  solution  makes  it 


possible,  by  employing  the  reciprocity  theorem,  to  derive  formulas 
for  calculating  the  radiation  characteristics  of  arbitrarily  ori¬ 
ented  electric  and  magnetic  dipoles  located  near  the  surface  of 
the  body. 


j 


. 


PART  3 


RECORDING  ELECTROMAGNETIC  RADIATION 


POSSIBLE  USE  OP  THE  MAGNETIC  MOMENT  OF  SEMICONDUCTOR 
PLASMA  ELECTRONS  FOR  RECORDING  CARBON  DIOXIDE  LASER  RADIATION 

G.  D.  Lobov,  V.  V.  Shtykov,  and  Ye.  I.  Gratsianskaya 

The  presence  of  a  nonlinear  effect  in  a  semiconductor  caused 
by  the  magnetic  moment  of  conduction  zone  electrons  was  indicated 
in  [48].  Also  made  in  this  publication  were  estimates  of  the 
magnitude  of  this  effect  in  the  millimeter  wave  band.  In  prin¬ 
ciple,  a  frequency  converter  based  on  this  effect  in  the  infrared 
band  can  be  constructed  using  semiconductors  with  small  effective 
mass.  From  the  results  obtained  in  [48]  it  follows  that  to  inten¬ 
sify  this  nonlinear  effect,  the  product  of  the  frequency  of  the 
effective  radiation  and  the  pulse  relaxation  time  coTq  of  the  semi¬ 
conductor  must  be  larger  than  unity.  This  condition  can  be  met  in 
the  infrared  band  at  room  temperature;  evidently  this  makes  it 
possible  to  solve  the  problem  of  the  uncooled  low-inertia  detector 
in  this  band. 

In  applications,  directing  the  permanent  magnetic  field  parallel 
to  the  propagation  of  electromagnetic  waves  may  prove  to  be  more 
convenient.  So  in  this  paper,  in  contrast  to  [48],  magnetization 
will  be  found  in  the  longitudinal  propagation  of  plane  electromag¬ 
netic  waves  in  a  semiconductor  plasma.  Since  in  actual  frequency 
converters  the  semiconductor  has  finite  dimensions,  when  we  set 
out  to  find  magnetization  we  take  into  account  reflections  from 
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the  interface.  Let  the  semiconductor  occupy  a  half-space . that 
has  the  boundary  z=0,  and  let  electromagnetic  waves  propagate 
along  the  z-axis.  Then  the  solution  of  the  differential  equa¬ 
tion  of  magnetization  presented  in  [48],  with  reference  to  the 
transmission  coefficient  T  of  the  electromagnetic  waves  in  the 
semiconductor,  will  be  of  the  form 

( 1 } 

are  the  constants  of  propagation  at 
the  corresponding  frequencies 
is  the  plasma  frequency  of  the  semi¬ 
conductor 

is  the  relaxation  time  of  a  pulse  in 
the  semiconductor 

is  the  effective  mass  of  the  conduc¬ 
tion  2one  electron  of  the  semiconductor 
is  the  dielectric  constant  of  the 
semiconductor  lattice 

is  the  index  of  refraction  of, the  semi¬ 
conductor  at  the  given  frequency 
is  the  attenuation  coefficient  at  the 
given  frequency 

In  Eq.  (1)  we  take  into  account  the  contribution  to  the  mag¬ 
netization  of  the  semiconductor  of  only  the  wave  for  which  the 
direction  of  rotation  of  the  electric  field  intensity  vector  coin¬ 
cides  with  the  direction  of  rotation  of  electrons  in  the  semicon¬ 
ductor. 


Using  the  results  obtained  In  [42],  the  square  of  the  modulus 
of  the  transmission  coefficient  can  be  represented  in  the  form 


ltf’«  [<><%•(* 


t:  : 


79 


4*  OS  <2  (6  10  H  4S 


where 


!•<* 


,  v  3*  ,  . 


is  the  cyclotron  frequency  of  the  semicon¬ 
ductor 


With  respect  to  Eq.  (2),  in  Eq.  (1)  we  isolate  the  part  dependent 

2  2 

on  the  concentration  of  the  conduction  zone  electrons  Uq/u  ,  the 
pulse  relaxation  time  u>Tq,  and  the  magnetic  field  intensity  (oj^/uj) 


S'), 


(3) 


where 


*Y 

ES 

ro1 


'/■-  //  azsrrviZTUZjCTTTSVS 


Presented  in  Figs.  1  and  2  are  the  above-indicated  functions. 

From  Fig.  1  we  see  that  the  maximum  value  of  M2(cu2)  rapidly  rises 

when  ojxn>3.  From  Fig.  2  can  be  found  the  value  of  the  permanent 
u  2  2 

magnetic  field  for  certain  u>tq  and  wQ/a)  .  From  a  comparison  of 

Figs.  1  and  2,  it  follows  that  to  obtain  the  maximum  value  of 

M2(o>2)  when  u)tq>3,  the  required  magnetic  field  is  approximately 

equal  to  its  value  at  the  point  of  cyclotron  resonance. 


The  variable  part  of  the  magnetization  of  the  semiconductor 
can  be  recorded  with  an  induction  coil  enclosing  the  semiconductor 
and  located  perpendicular  to  the  external  magnetic  field.  Let  us 
estimate  the  parameters  of  the  frequency  converter  at  A=10.6  micro¬ 
meters,  in  which  a  semiinfinite,  longitudinally  magnetised  plasma 
cylinder  with  radius  R>>X  is  used. 

The  above-named  condition  allows  us  to  use  the  results  derived 
above  in  making  estimates  of  the  converter  parameters.  If  the  re¬ 
cording  induction  coil  directly  surrounds  the  semiconductor  cylin¬ 
der,  the  open-circuit  voltage  at  coil  ends  is  [50] 

I'll  (1° 
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where  £— f  is  the  cross  sectional  area  of  a  coil  turn  and  ALm 
is  the  number  of  coil  turns.  If  the  laser  beam  illuminates  the 
entire  plasma  cylinder  cross  section,  then  Eq.  (4),  with  reference 
to  Eq.  (3)  takes  on  the  form 

(5) 

In  [5]  it  is  considered  that  the  intensity  of  electromagnetic 
waves  incident  at  the  interface  is  ,  where  Wm'{%  is  the 

characteristic  impedance  in  vacuo,  and  S  is  the  beam  cross  sectional 
area. 


Estimate  calculations  showed  that  at  the  output  of  the  fre¬ 
quency  converter  at  A=10.6  micrometers  (C0o  laser),  for  InSb  with 

—  ^  2  ^ 

m*=l. 52*10  kg,  0^=16,  <dt0=10,  the  latter  corresponds  to 

U=0.6  m2/V»s  [51],  when  P1=P2=1  W  and  to  the  coil  turn  number 

M  _  =10;  the  open-circuit  voltage  is  approximately  5  microvolts 

and  Is  attained  for  a  magnetic  field  Intensity  of  160,000  gauss. 

To  lower  the  magnetic  field  Intensity,  semiconductors,  with  a  lower 

effective  mass  can  be  used.  If  an  HgTe  semiconductor  or  a  ternary 

compound  CdxHg^_xTe  is  used  as  the  working  material  In  the  above- 

described  converter,  with  an  effective  conduction  zone  electron 

mass  of  0.003m  ,  the  open-circuit  voltage  for  these  parameters  will 

“  P 

be  approximately  26  microvolts,  for  a  mobility  y=2.5  m  /V-sec,  and 
a  magnetic  field  intensity  of  31.6  kilogauss. 

Thus,  utilizing  cdxHSi_xTe  in  this  converter  is  preferred  over 
InSb.  Synthesizing  semiconductor  compounds  with  an  even  lower  ef¬ 
fective  mass  makes  it  possible  to  substantially  lower  the  required 
magnetic  field  and  to  build  detectors  for  \=10.06  micrometers  with 
magnetic  fields  realistically  attainable  at  the  present  time  without 
cooling. 
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CONVERSION  OF  C02  LASER  RADIATION  TO  MILLIMETER  BAND 

RADIATION  USING  THE  MAGNETIC  MOMENT  OF  CONDUCTION  ZONE 
ELECTRONS  OF  SEMICONDUCTORS 

V.  V.  Shtykov 

The  problem  of  the  conversion  of  infrared  radiation 
frequency  in  a  semiconductor  placed  In  a  constant  magnetic  field 
Is  examined  in  [113].  An  induction  coil  Is  suggested  In  recor¬ 
ding  the  variable  magnetic  moment,  varying  with  time  with  the 
difference  frequency  co^.  But  there  is  some  interest  in  Investi¬ 
gating  phenomena  when  infrared  radiation  is  converted  to  milli¬ 
meter  band  radiation,  since  the  frequency  can  be  1010  to  1011 
hertz  (see  equation  for  the  magnetization  vector  M  in  [48]).  The 
high  value  of  the  difference  frequency  precludes  the  use  of  an 
Induction  coil  in  recording  the  variable  part  of  M.  In  this  case 
a  system  suitable  for  microwave  applications  must  be  used;  in 
the  system  an  electromagnetic  wave  with  the  frequency  will  be 
excited. 

Let  us  examine,  by  way  of  example,  a  variant  of  an  infrared 
radiation  frequency  converter  in  which  a  rectangular  metal  wave¬ 
guide  is  used  (see  Fig.  1). 

To  solve  the  problem  of  wave  excitation  In  a  waveguide  we  will 
assume  that  the  nonlinearity  in  the  semiconductor  is  small  [48]. 
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This  allows  us  to  regard  the  magnitude  of  the  variable  part  of 
the  vector  M  as  a  given  function  and  set  it  in  correspondence 
to  the  outside  magnetic  flux  with  the  complex  amplitude 

j  (1) 

The  complex  amplitude  of  the  variable  part  (u>£  ,  which  appears 
in  Sq.  (1),  was  found  for  two  limiting  cases  of  the  propagation 
of  electromagnetic  waves  with  frequencies  and  with  respect 
to  the  vector  of  the  constant  magnetic  field  HQ  in  [48]  and  [113]. 

Prom  the  results  of  these  studies  it  follows  that  if  the 
principal  waves  are  traveling  waves,  then  is  also  a  traveling 

wave.  This  circumstance  imposes  a  number  of  features  on  the  con¬ 
version  of  infrared  radiation  to  a  millimeter  band  wave.  It  turns 
out,  for  example,  that  when  several  phase  relationships  are  satis¬ 
fied,  waves  from  individual  regions  of  the  semiconductor  may  be 
added  together;  this  may  lead  to  a  larger  conversion  coefficient. 
This  phenomenon  is  an  analog  of  phenomena  that  are  discussed  in 
studies  on  nonlinear  optics  [54]. 

If  we  consider  that  the  component  A7t(c*>j)  parallel  to  the  field 
Hq  has  a  maximum  magnitude,  and  that  for  the  most  efficient  wave 
excitation  the  vector  5"  must  be  parallel  to  the  vector  H  of  the 
excited  wave  [10],  Hq  must  be  oriented  parallel  to  H.  When  the 
principal  type  of  oscillations  in  a  rectangular  waveguide  is  used, 
the  last-named  condition  can  be  satisfied  by  orienting  HQ  either 
perpendicular  to  the  narrow  wall  of  the  waveguide  or  along  its 
longitudinal  axis.  If  also,  as  In  [113] ,  we  consider  the  case  of 
longitudinal  propagation  of  the  principal  waves,  the  variant  shown 
in  Fig.  1  makes  it  possible  to  use  longer  waveguide  specimens. 

By  employing  the  method  outlined  in  [10]  and  Eq.  (1),  we  can 
find  the  complex  amplitude  of  the  electric  field  at  the  waveguide 
center  in  the  cross  section  z=L 

Em  (**>»)  ***  Y^K)***- ^ 


(2) 


In  deriving  Eq.  (2)  we  neglected  the  effect  of  the  semiconductor 
wafer  for  waves  in  the  waveguide,  since  a>>A.  In  the  expression 
derived  the  attenuation  of  waves  in  the  waveguide  was  allowed  for 
by  introducing  the  complex  longitudinal  wave  number  h=h'-jh". 

If  waves  with  circular  polarization  propagate  in  the  semicon¬ 
ductor,  the  direction  of  rotation  of  vector  E  of  the  waves  coin¬ 
ciding  with  the  direction  of  electron  spin,  then  /fyfu/,)  can  be  re¬ 
presented  in  the  form 

Assuming  that  Mzg  does  not  depend  on  x  and  y,  by  performing  inte¬ 
gration  in  Eq.  (2)  with  respect  to  z  from  0  to  L,  we  get 


The  dependence  of  the  field  ExQ  on  the  length  of  the  wave  interac¬ 
tion  region  can  be  described  with  a  function  conveniently  repre¬ 
sented  in  the  complex  plane  (see  Fig.  1),  if  we  transform  it  to 
the  form 

f(Q  ■  «* pC’aJC*)-  exp  ~  (3) 

by  Introducing  the  parameter  -h')l»AKL  . 

From  Fig.  1  it  follows  that  when  ««♦*,<  ah  ,  f(L)  is  oscil¬ 
latory,  and  when  o t,*«Lt>xK  ,  the  oscillations  are  weakly  manifested 
Emphasis  is  given  to  the  case  *  In  this  case,  f(L)  is  a 

real  variable.  This  regime  corresponds  to  the  condition  of  phase 
matching,  discussed  in  studies  on  nonlinear  optics. 

From  Fig.  1  it  follows  that  the  deviation  from  this  regime  is 
not  large  when  404*  .  When  this  condition  is  satisfied, 

f(L)*»xp  (-KL)  -  exp  C-  «  ♦«*)  U  (4) 

Eq.  ( 4 )  attains  a  maximum  value  when 

where 


The  field  amplitude  Ex0  in  this  case  attains  a  maximum  value,  equal 


to 


where 


Shown  in  Fig.  2  are  h‘LM  and  F(Z)  as  functions  of  the  para¬ 
meter  Z.  From  the  figure  it  follows  that  it  is  best  to  function 
when  sTTcTt —  *°.i  1  which  corresponds  to  the  length  Lm<  $  . 

9 

Corresponding  to  the  value  Exq  found  when  *#*afi*<*  is  the  in¬ 
tensity  of  waves  with  a  difference  frequency  equal  to 

9 

Using  the  expression  from  Mz0  obtained  in  [113]  and  considering 
2*n.  ,  where  a+-  is  an  element 

of  the  conductivity  tensor;  when  h’>>h",  the  expression  for 
P(u>3;Lra)  can  be  represented  in  the  form 

PH*  A.)  •  p< p-  F‘m  <  s ) 

Here  and  P2  are  the  intensities  of  the  principal  waves  in  the 
crystal,  Bg  is  electron  mobility,  and  Wg  is  the  characteristic 
impedance  of  the  waveguide. 


Let  us  estimate  the  intensity  of  difference-frequency  waves 
when  radiation  from  a  C02  laser  (  )  is  converted  to  milli¬ 

meter  band  radiation  with  wavelength  equal  to  approximately  5.6 
mm  (see  [115]).  When  a  standard  4  millimeter  waveguide  is  used 
and  with  Z<0.01,  the  Intensity  PCw^)  is 

Pfaipm)  *  PiP* 

If  an  InSb  semiconductor  crystal  is  used  as  the  material,  77 
2 

Be  can  reach  70  m  /s  at  the  temperature  of  liquid  nitrogen.  When 

P.=P-=100  W,  which  is  quite  realistic  for  a  C0?  laser,  the  aiffer- 

1  d  -10  £ 
ence-frequency  intensity  will  be  5 ‘10  W. 
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Thus,  there  is  the  opportunity  of  the  experimental  observa¬ 
tion  of  the  conversion  of  coherent  infrared  radiation  to  milli¬ 
meter  band  radiation  using  the  magnetic  moment  of  the  semicon¬ 
ductor  conduction  zone  electrons. 


CONVERTER  OP  FREQUENCY  OP  COHERENT  INFRARED  RADIATION  TO  THE 
MICROWAVE  BAND  USING  A  WAVEGUIDE  WITH  AN  ELECTROOPTICAL  CRYSTAL 


V.  I.  Bogatkin 

The  conversion  of  the  frequency  of  coherent  infrared  radia¬ 
tion  to  the  microwave  band  can  be  attained  by  the  excitation  of 
a  waveguide  filled  with  an  electrooptical  crystal  with  nonlinearly 
polarized  current  at  the  difference  (intermediate)  frequency 
cJm •  u)r - o»c  Induced  when  heterodyne  waves  and  Infrared  band  sig¬ 
nals  are  propagated  in  the  crystal  with  the  frequencies  Wp  and 
a)  ,  'respectively  [52,56]. 

Below  is  outlined  a  method  of  calculating  the  parameters  of 
these  converters  with  reference  to  the  structure  of  the  fields  of 
all  three  interacting  waves  and  losses  in  the  crystal;  additionally, 
the  transmission  coefficient  is  calculated  with  respect  to  the  ex¬ 
periment  presented  in  [52].  The  main  symbols  correspond  to  those 
adopted  in  [56]. 


The  complex  amplitude  of  the  electric  field  of  the  interme¬ 
diate-frequency  wave  in  the  waveguide  can  be  described  in  the 


form  [53] 


Here  the  function  t^(x, y)  describes  the  structure  of  the  field  of 
the  excited  wave  type  in  the  waveguide  cross  section,  and  the 
3lowly  varying  amplitude  (a)  obeys  the  equation 
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(1) 


SjW„  • 

S  • 

where  _  . 

^■i^iuj:u,d.dy 

is  the  normalized  energy  flow  in  the  waveguide  (neglecting  the 
spatial  dispersion  of  the  crystal,  see  [53]).  Integration  is 
carried  out  over  the  waveguide  cross  section.  The  coefficient 
<*„  is  the  attenuation  coefficient  with  respect  to  the  field  of 
the  excited  wave  with  allowance  for  losses  in  the  crystals  and 
in  the  waveguide  walls.  The  tensor  describes  the  density  of 
the  nonlinearly  polarized  current  induced  in  the  crystal  when 
acted  on  by  the  fields  of  the  heterodyne  wave  and  of  the  signal 
wave  with  the  amplitudes 

K and  Ec- 

respectively  [53,54,551.  In  Eq.  (1),  A**xr-*f-*4  .  Equations 
analogous  to  (1)  obtain  also  for  the  amplitudes  of  the  heterodyne 
wave  and  of  the  signal  wave;  when  combined  with  Eq.  (1),  they 
describe  the  Interaction  of  waves  with  three  frequencies. 

For  a  small  transmission  coefficient  of  the  converter  with 
respect  to  power,  Eq.  (1)  can  be  solved  by  assuming  the  quantities 

9  • 

Ur  and  U  to  be  assigned.  Considering  the  losses  in  the  crystal, 

we  can  write  ,  where  the  quantities  Up  and 

U  apply  to  some  initial  crystal  cross  section.  In  this  cas§  the 
G0 

solution  of  Eq.  (1)  coincides  with  the  solution  obtained  from  Eqs. 
(1)  and  (6)  in  [5],  if  the  quantity  mf„  in  [5]  is  replaced  with 

(U>J d* d*  (3) 

and  if  .  In  particular,  all  results  applying  to  the  optimal 

crystal  length  and  the  transmission  coefficient  as  a  function  of 
the  parameters  (introduced  in  [56])  derived  from  analyzing  the 
idealized  case  of  the  interaction  of  plane  waves  are  valid  for 
conversion  in  the  waveguide.  The  results  of  the  width  of  the 
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radiation  pattern  in  the  passband,  obtained  in  [56],  can  also  be 
extended  to  the  case  of  conversion  in  the  waveguide. 

Now  let  us  calculate  the  transmission  coefficient  of  the  con¬ 
verter  as  applied  to  the  experiment  conducted  in  [52].  In  it,  the 
radiation  from  a  CC^  laser,  containing  oscillations  at  two  frequen¬ 
cies,  was  focussed  with  a  long-focal  length  lens  at  the  center  of 
a  rectangular  waveguide  filled  with  GaAs  crystal  and  excited  H10 
waves.  By  setting  the  waveguide  cross-sectional  dimensions  equal 
to  a  and  b,  with  respect  to  the  x-  and  y-axes,  respectively,  we  can. 
describe  the  fields  of  laser  radiation  in  the  crystal  in  the  form 


ir  ■  if.,-'*  a  eV.,. 


where  R  is  the  effective  radius  of  the  laser  beam  entering  the 
crystal  in  the  z=0  cross  section.  Strictly  speaking,  these  func¬ 
tions  cannot  describe  the  field  of  waves  propagating  in  the  wave¬ 
guide,  but  when  R<<a,b  and  assuming  that  the  geometrical  optics  is 
valid  for  the  laser  radiation  in  the  waveguide,  they  are  a  good 
approximation.  The  eigenfunctions  for  the  excited  H1Q  wave  we  can 
write  in  the  form  .  — f.  3* 


where  is  the  characteristic  impedance  of  the  waveguide  filled 
with  an  isotropic  GaAs  crystal.  To  construct  Eq.  (1),  let  us  write 
/T.  according  to  Eq.  (2) 


Let  us  compute  the  integral  in  Eq.  (3) 


where 


m  .  0  /r  iY 
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(3) 


•When  R<<a,b,  Eq.  (6)  can  be  easily  computed: 

jfCHC « 

By  multiplying  Eq.  (1)  by  Eg  and  inserting  in  it  Eqs.  (5)  and 
(8),  we  get  an  equation  for  the  complex  amplitude  of  the  field 
E  of  the  H1Q  wave  in  the  waveguide  center 


(9) 


If  by  selecting  the  waveguide  dimension  a  the  condition  Ak=0  is 
satisfied  and  if  the  crystal  length  is  small,  such  that  we  can 
neglect  the  attenuation  of  all  waves,  then  by  assuming  in  Eq.  (9) 
that  ctrmmka given  the  boundary  condition  a,-- 0  we  get  the 
solution 

£(a)  v  40 

do) 


Considering  that  the  power  £(*)  transported  by  the  excited  wave 
in  the  z  cross  section  is  given  by  Eq.  (5),  when  in  place  of 
Eg  we  must  substitute  E(z),  and  also  that  the  amplitudes  Ep 

and  E  are  associated  with  the  powers  Pr  and  P  of  the  hetero- 
0  0  c0 
dyne  wave  and  the  signal  wave  within  the  crystal  at  the  converter 

input  by  the  relationships 


£ 


l£  r‘ 


where  Wp ,  W  are  the  characteristic  impedances  of  the  crystal  for 
the  corresponding  waves,  we  get  the  transmission  coefficient  with 
respect  to  power  when  the  crystal  length  is  z: 


4®. 


(11) 


Let  us  compute  the  nonlinearity  parameter  m  in  Eq.  (11). 

The  components  of  the  nonlinear  polarization  vector  for  the  GaAs 
crystal  in  the  crystallographic  coordinate  system  x' ,y’ ,z'  are 
equal  to  [54,55] 
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In  conclusion  we  present  plots  (Pigs.  1-4)  characterizing 
the  change  in  power  of  the  intermediate-frequency  wave  at  the 
converter  output  when  the  crystal  length  changes,  with  allowance 
for  losses  and  the  mismatch  of  phase  velocities.  The  plots  were 
constructed  on  the  basis  of  the  formulas  derived  in  Eq.  (5)  when 
an  analysis  was  made  of  the  interaction  of  plane  waves.  They 
were  normalized  along  the  x-axis  to  the  optimal  crystal  length, 
and  along  the  y-axis,  to  the  transmission  coefficient,  given  an 
optimal  crystal  length.  The  symbols  correspond  to  those  adopted 
in  [56]. 


CHARACTERISTICS  OP  THE  CONVERSION  OP  COHERENT  INFRARED 
RADIATION  TO  THE  MICROWAVE  BAND  IN  A  RESONATOR  CONTAINING 
AN  ELECTROOPTICAL  CRYSTAL 

V.  I.  Bogatkin 

Conversion  of  the  frequency  of  coherent  infrared  radiation 
to  the  microwave  band  can  be  attained  when  an  electrooptical 
crystal  is  placed  not  only  in  a  waveguide  [52,56,115],  but  only 
in  a  microwave  resonator  [56].  Certain  properties  of  these  con¬ 
verters  are  Investigated  In  this  article  by  analysis  of  Ideal¬ 
ized  cases.  The  symbols  and  the  subscripts  correspond  to  the 
symbols  and  subscripts  adopted  In  [115]. 

The  power  transmitted  by  nonlinearly  polarized  currents 
with  the  frequency  «»)„  to  the  excited  type  of  resonator  oscilla¬ 
tions  is 

dV-  (i) 

here  integration  Is  carried  out  with  respect  to  the  resonator 
volume , 

(2) 

is  the  density  of  the  nonlinearly  polarized  current; 

is  the  complex  amplitude  of  the  electric  field  vector  of  the 

excited  oscillations;  is  the  function  describing  the 


the  spatial  distribution  of  the  field.  The  amplitude  cofactor 
V*  ,  assuming  that  equals  the  resonance  frequency  of  the 
resonator,  is  given  by  the  expression  [53] 


where 


(3) 

(4) 


is  the  normalized  energy  of  the  resonator  field;  Qm  is  the  Q 
of  the  loaded  resonator. 


By  substituting  Eq.  (2)  in  Eq.  (3),  and  then  Eqs.  (2)  and 
(3)  in  Eq.  (1) ,  we  get 

•l-tifrUOK&f  (5) 


Eq.  (5)  yields  the  total  power  transported  by  the  nonlinearly 
polarized  currents.  Part  of  this  power  is  dissipated  in  the  re¬ 
sonator  owing  to  the  losses,  and  only  part  of  it  can  be  used  in 
the  form  of  a  signal  in  the  load  associated  with  the  resonator. 

We  can  show  that  the  maximum  power  released  in  the  load  is 


I 


where  Qq  is  the  Q  of  the  unloaded  resonator.  Therefore 


(6) 


Let  the  resonator  be  formed  of  two  infinite  planes  whose 
enclosed  space  is  filled  with  an  isotropic,  nonlinear  dielectric 
with  thickness  l  .  Let  us  assume  that  the  plane  TEM  wave  with 
field  components  £mm 


is  excited,  where  z  is  the  coordinate  measured  along  the  z— axis 
perpendicular  to  phe  resonator  walls;  *nt*xn  ,  n  -1,2,3...  ; 

the  unit  vectors  and  I  are  perpendicular.  Let  us  calculate 
the  quantity  in  .Eq.  (4): 

Wm-ts  E'cos'kj)  cfcm 

*  4  m  C*  t  ( 8 ) 

Here  S  is  the  cross-sectional  area  of  the  resonator.  Physically 
speaking,  it  is  reasonable  to  calculate  the  power  received  per 
unit  resonator  cross-sectional  area,  so  in  the  following  we  will 
assume -S^  to  be  equal  to  the  unit  area.  When  calculating 
Eq.  (8),  it  is  also  assumed  that  and  do  not  depend  on 
frequency  near  <*>/». 

Excitation  of  Currents  by  a  Traveling  Wave 

We  assume  that  plane  heterodyne  and  signal  waves  of  the  TEM 
mode  propagate  through  the  resonator  along  the  z-axis.  In  this 
case,  we  assume  the  resonator  walls  to  be  transparent  to  these 
-waves.  To  some  extent  this  case  corresponds  to  the  actual  situ¬ 
ation,  when  signal  and  heterodyne  beams  enter  and  exit  through 
holes  in  the  resonator  walls.  The  wave  amplitudes,  with  allowance 
for  attenuation,  we  write  in  the  form 

4.tt«**V**»  (9, 

By  inserting  Eqs.  (8)  and  (9)  in  Eq.  (6)  we  get  the  power  received 
per  unit  resonator  area: 

K  K  nAn,.  e> 


where  *'*)•"*«  •  are  characteristic 

impedances.  e»  are  the  Poynting  vectors  of  the  correspon¬ 

ding  waves  after  their  entrance  Into  the  resonator  (when  z=0). 

The  function  f  is  equal  to  the  following  expression: 

I/i  -  -  *■  n  n 


Lv  _  ^  -  H)  Ze  cot 
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Let  us  determine  the  conditions  under  which  Eq.  (11)  is  at  a  max¬ 
imum;  in  view  of  the  complexity  of  function  f,  let  us  limit  our¬ 
selves  to  an  approximate  analysis.  The  expression  in  the  square 
brackets  in  the  denominator  of  Eq.  (11)  has  a  minimum  with  change 
in  k,  which  is  attained  when 

(12) 

The  numerator  has  a  maximum  when  the  term  in  it  is  positive  and 
JcoswljM  .  This  requires  that  •?»  ,  where  q  is  a  positive 
integer.  By  equating  both  values  to  k,  we  find  that  we  must  have 

in  order  that  the  denominator  minimum  and  the  maximum  of  the 
numerator  in  Eq.  (11)  coincide. 


In  calculating  q,  it  was  assumed  that  Prom  the  latter 

it  will  be  clear  that  aH^l  is  necessary  in  order  to  maximize  f, 
therefore  q<<l. 


So,  given  the  condition  •/I*,,  f  ,  we  cannot  carry  out  Eq.  (12) 
and  q=l  at  the  same  time.  So  when  k  changes,  the  maximum  of  the 
numerator  and  the  minimum  of  the  denominator  of  Eq.  (11)  do  not 
coincide.  Further  analysis  shows  that  under  condition  (12)  the 
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value  of  f  exceeds  by  about  one  order  of  magnitude  its  value 
when  .  So  we  assume  that  Eq.  (11)  is  at  a  maximum 

when  Km*n  j  then  we  get 

m  ifct  (*-**)'  (13) 

A  plot  of  function  (13)  is  given  in  Fig.  1  (curve  1).  The  curve 
has  a  maximum  when 

oil  •  4,2S , 


(14) 


is  approximately  equal  to  0.1. 


Thus,  the  maximum  possible  transmission  coefficient  of  the 
converter  under  study  is 


(15) 


Excitation  of  Currents  with  a  Standing  Wave 


Let  us  assume  that  in  the  resonator  under  study  the  hetero¬ 
dyne  and  signal  waves  are  reflected  from  the  internal  wall  when 
z-l  and,  on  propagating  in  the  opposite  direction,  exit  from 
the  resonator  when  z=0.  Then  for  the  heterodyne  and  signal  waves 
we  will  have  the  expression 

The  minus  sign  in  front  of  the  second  terms  takes  into  account 
the  loss  of  a  half  wave  at  reflection. 

When  computing  the  components  in  Eq.  (6) ,  we  will 

consider  only  the  products  of  the  first  and  second  terms  in 
Eqs.  (16)  and  (17).  The  rejected  terms  correspond  to  the  weak 
interaction  of  the  heterodyne  and  signal  waves  moving  in  opposite 
directions.  By  calculating  Eq.  (6),  we  arrive  at  an  expression 


(16) 

(17) 


analogous  to  Eq.  (10): 


where 

•  [<+ e“*i  4r«’*^ 

• 

Analyzing  this  function  at  the  maximum,  like  the  above  analysis, 
again  leads  to  the  condition  of  phase  matching  *»#«.  Here  the 
function  f,  is  transformed  to 

(20) 

Its  plot  is  given  in  Pig.  1  (curve  2).  It  attains  a  maximum 
equal  to  0.072  when  a.l-2.5.  The  maximum  transmission  coefficient 
in  this  case  is 

Sit**".  (n) 

By  comparing  Eq.  (21)  with  Eq.  (15)  we  discover  that  using 
the  reflected  heterodyne  and  signal  waves  led  to  a  reduction  of 
the  transmission  coefficient  by  a  factor  of  1.4  and  to  an  increase 
in  the  optimal  resonator  length  by  a  factor  of  2.  This  result 
is  explained  by  the  fact  that,  when  reflected  by  the  resonator 
wall,  both  waves  (heterodyne  and  signal)  change  phase  by  tt. 

Since  Eq.  (2)  includes  the  products  of  wave  amplitudes,  the  re¬ 
flected  waves  will  produce,  near  the  wall,  a  nonlinear  polariza¬ 
tion  vector  in  the  same  phase  as  the  incident  waves.  The  field 
of  the  excited  oscillation  (7)  is  formed  by  two  waves  moving  in 
opposite  directions  and  varying,  at  reflection  by  the  wall,  in 
phase  by  it.  As  a  result,  if,  for  example,  the  heterodyne  and 
signal  waves  moving  in  the  positive  direction  of  the  z-axis  give 
up  their  energy  to  the  intermediate-frequency  wave  moving  in  the 
same  direction,  after  reflection  of  all  three  waves  the  energy 
will  be  pumped  in  the  opposite  direction  because  of  the  disturbed 


(18) 


(19) 
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matching  of  the  phases  of  nonlinearly  polarized  currents  and  the 
electric  field.  So  the  conversion  efficiency  suffers.  The  pre¬ 
sence  of  losses  for  the  heterodyne  and  signal  waves  compensates 
the  effect  of  the  reverse  transmission  of  energy,  leading  to  an 
increase  in  the  optimal  resonator  length. 


Let  us  assume  that  the  effect  of  phase  discontinuity  at  re¬ 
flection  is  eliminated.  In  actual  conditions  this  can  be  achieved 
by  making  a  hole  in  the  resonator  wall  at  z=&  and  positioning  the 
mirror  reflecting  the  heterodyne  and  signal  beams  in  a  location 
outside  the  resonator  at  a  distance  of  one-quarter  intermediate- 
frequency  wavelength  in  the  dielectric.  When  making  the  calcula¬ 
tions  for  this  case,  a  plus  sign  must  be  placed  at  the  second 
term,  for  example,  in  Eq.  (17)  (we  neglect  the  attenuation  in  the 
dielectric  section  projecting  from  the  resonator).  Then  in  Eq. 
(18),  in  place  of  f . ,  the  function 


(22) 


will  appear.  As  before,  it  will  have  a  maximum  at  equal  to 


f% 9 2 is  ft"***/*  zf(t*Q 


(23) 


By  comparing  Eqs.  (23)  and  (13)  we  see  that  when  there  is 
reflection  with  compensation  of  the  half-wave  loss,  the  optimal 
resonator  length  is  reduced,  and  the  transmission  coefficient 
is  increased  by  a  factor  of  2  compared  with  the  case  of  excita¬ 
tion  with  a  traveling  wave. 


We  note  that  compensation  of  the  half-wave  loss  is  not  re¬ 
quired  when  there  is  excitation  by  TM  waves  (for  example,  in  a 
cylindrical  resonator  [56]). 


Comparison  of  Transmission  Coefficients  in  Waveguide  and 
Resonator  Converters 


If  losses  in  the  resonator  are  caused  only  by  losses  in  the 
dielectric,  we  can  find  that 

w,- 

where  is  the  attenuation  coefficient  of  the  amplitude  of  the 

It 

plane  wave  with  frequency  cO  in  the  dielectric. 

n 


The  maximum  transmission  coefficient  in  the  case  of  excita¬ 
tion  with  a  standing  wave  with  compensation  for  half-wave  loss 
will  then  be  written  as 


(25) 


We  can  show  that  the  idealized  waveguide  converter  (examined  in 
[56])  with  a  plane  wave  has  a  transmission  coefficient  equal  to 


(26) 


where  $(*,**')  is  a  function  in  Eq.  (12),  in  [56]. 


Let  us  inspect  the  ratio 

/ 

Mmf '(*,&*')  (27) 

where  •  When  »•*  and  Ak'=0,  the  gain  (26)  in  the 

transmission  coefficient,  for  a  resonator,  is  approximately 
1.5.  Calculations  show  that  the  optimal  resonator  length  in 
this  case  is  also  shorter  than  the  waveguide,  by  a  factor  of  1.5. 


From  Eq.  (27)  it  follows  that  GJJ  can  exceed  GMb  without 
bound,  if  ae—  0. 
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RECORDING  SUBMILLIMETER  AND  INFRARED  RADIATION  WITH  THE 
TUNNEL  METAL-TO-DIELECTRIC-TO-METAL  (MDM)  JUNCTION 

G.  D.  Lobov  and  A.  N.  Nenashev 

Characteristics  of  Tuhnel  Current  for  Weak  Signals 

Current  flowing  by  tunnel  penetration  through  the  potential 
barrier  in  the  MDM  contact  is  defined  as  follows  [57,58]: 


here  e  and  m  are  the  electron  charge  and  mass;  k  is  the  Boltzmann 
constant ;n  is  the  Fermi  energy;  T  is  temperature  in  degrees  Kelvin; 
v  is  the  applied  voltage;  h  is  the  Planck  constant;  E  is  the 

A 

fraction  of  the  electron  energy  in  the  direction  perpendicular 
to  the  potential  barrier;  Em  is  the  maximum  electron  energy; 

S  Is  the  contact  area;  l  is  the  dielectric  thickness; 

2(u.mpC-/ifm»‘ru*c8j  Is  the  probability  of  tunnel  pene¬ 
tration  of  an  electron  with  energy  Ex;  $(x)  is  the  shape  of  the 
potential  barrier;  and 

When  T>0K,  the  representation  of  D(E  )  in  the  form  [59] 
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If  1-c*T>(2f)\T  ,  we  can  neglect  the  third  term  in  Eq.  (2). 

In  [58]  it  is  shown  that  this  condition  is  usually  satisfied 
and  the  coefficients  b  and  c  are  defined  by  the  following  expres¬ 
sions  : 

fatty)  (3) 

If  the  voltage  iS  Vutinuit  is  applied  to  an  MDM  system,  the  detected 
signal  is  determined  by  the  following  expression  [57]: 

(4> 

where  R_y  is  the  contact  impedance  for  the  tunnel  current. 

Using  the  representation  of  D(E  )  in  the  form  of  Eq.  (2) 
and  expanding  Eq.  (1)  in  a  series  in  terms  of  v,  for  the  coeffi¬ 
cients  in  Eq.  (4)  we  get  the  expressions 


/  mJr  {  £  £mrM)(*«Tc-) 

♦  '  (x*7c^ cosecXxxTc^jDl-unfxHTc^lr 


Cm-C 


(xkTcJ  ciy  (ttTqJ)}, 

R  -  *>c"  HsJSSSa 

•  sr *Tcm 


(5) 

(6) 

(7) 


In  Eqs.  (5)-(7),  the  first  digits  in  the  subscripts  at  the  coef¬ 
ficients  b  and  c  correspond  to  the  following:  "2"  is  the  forward 
characteristic  and  "l”  is  the  reverse  characteristic.  The  second 
digits  give  the  order  of  the  derivative  of  the  expressions  in 
Eq.  (3)  in  terms  of  the  argument  ev. 


Allowing  for  the  Effect  of  Mirror  Image  Forces 

We  will  adopt  the  expression  for  the  potential  barrier  with 
allowance  for  the  mirror  image  forces,  as  in  [60],  in  the  form 
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(3) 


<P(*)-<PT(x)  *Vi. 

where  K:  *  *  J7yf)~  5  A*  \Jfgf  5  £  is  the  relative  dielectric  con¬ 
stant;  and  $T(x)  is  the  shape  of  the  trapezoidal  barrier. 


Since  it  is  difficult  to  derive  expressions  for  b  and  c 

in  final  form  from  Eq.  (3)  with  reference  to  Eq.  (8),  let  us 

derive  an  approximate  expression,  considering  that  X  is  a  small 

quantity.  Expanding  Eq.  (3)  in  series  in  terms  of  X  ,  after 

substituting  Eq.  (8)  in  it  and  neglecting  the  terms  proportional 
2 

to  x  ,  we  get 


($«  ♦  #«)  * 

>■  A  .  &  ] 

*  A  ^  a  mx  A  I  K.  rA71  1 


(9) 


•1  li  T  T, 

♦  j 


Here 

w 1 


The  coefficients  for  the  corresponding  trapezoidal  barrier  are 
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•  a  (<:  -  u  <v.  <  ■  -a  (c ;-u  *;*), 


(10) 


Fig.  1. 

(1)  Contact  temperature  0°K 

(2)  Contact  temperature  300°K 

Shewn  In  Fig.  1  (1)  1?  the  rectified  voltage  as  a  function  of 
the  dielectric  thickness,  derived  from  Eqs.  (4)— (10) .  From  the 
functions  shown  It  Is  clear  that  the  effect  of  the  mirror  image 
forces  is  quite  large  and  even  when  e<8,  a  signal  with  either 
positive  or  negative  polarity  can  be  obtained. 

Allowing  for  the  Thermionic  Emission  Current 

If  the  contact  temperature  of  an  MDM  is  higher  than  0°K, 
within  the  metal  there  will  be  electrons  whose  energy  is  larger 
than  the  work  yield  function  <$.  On  leaving  the  metal,  these 
electrons  produce  the  thermionic  emission  current  [ 6 1 ] 


In  .i*3rh‘%m»K'TlS[exp  f*f)  -  exp(--f)]  • 

•I,  C**P  (-%?)  -  ex p  (-£?)], 

where  and  <p^m  are  the  potential  barrier  maxima  in  metals  over 
the  Fermi  levels. 


If  the  voltage  v  is  applied  between  the  electrons,  then 

9im  "  9m  *  ^ 

and 

(’&)£<- **p  (‘  *r)  J  ^11) 

By  defining  4>lm  analogous  to  [61]  for  an  electrode  with  a  smaller 
work  function,  negatively  charged,  we  find  that  [61] 

^  ♦  efop-toPe'f*]  ( ]_2 ) 

If  the  electrode  with  a  larger  work  function  is  negatively  charged, 
then 

9Hmm  {i*ri**f{(*9*ei»] ff*]*  (13) 


By  expanding  Eq.  (11)  in  series  in  v,  using  Eqs.  (12)  and  (13), 
we  will  have  the  following  expression  for  the  rectified  voltage 
due  to  thermionic  emission 


(14) 


where 


& 


~a3«*n 


C*C(e* 


(15) 


By  representing  and  as  current  generators  with  internal 
impedance  fir*  and  Bn  ,  for  the  total  voltage  at  the  contact 
s  )  we  get  the  expression 


(16) 


Shown  in  Pig.  1  (2)  are  the  results  of  calculation  based  on 
Eq.  (16)  with  the  involvement  of  Eqs.  (4),  (7),  (14),  and  (15) 
for  the  contact  temperature  T=300°K.  Prom  the  figure  it  is  clear 
that  up  to  a  certain  dielectric  thickness  the  effect  of  temperature 
is  small.  But  when  the  dielectric  thickness  is  increased,  corres¬ 
ponding  to  an  increase  in  contact  impedance,  a  second  change  in 
the  polarity  of  the  rectified  voltage  is  possible. 

Experimental  Results 

The  experiments  were  conducted  with  a  point  contact  between 
a  15  micrometer  tungsten  wire  and  an  aluminum  base  coated  with 
a  natural  oxide  layer.  The  tungsten  wire  was  sharpened  electro- 
lytically  and  the  tip  was  monitored  with  a  microscope.  The  mech¬ 
anical  system  made  it  possible  to  finely  regulate  the  wire  pressure 
on  the  base,  by  varying  the  thickness  of  the  dielectric  film.  The 
relative  changes  in  the  dielectric  thickness  were  monitored  by  the 
change  in  the  contact  at  direct  current.  The  following  device  was 
used  as  supplies  of  microwave  oscillations:  a  klystron  operating 
in  the  5  mm  range,  and  a  CC^  laser.  In  both  cases  amplitude  modu¬ 
lation  of  the  radiation  was  employed.  In  the  experimental  part  of 
the  study,  the  magnitude  and  shape  of  detected  pulses  were  investi¬ 
gated  as  a  function  of  dielectric  film  thickne'ss. 

Shown  In  Fig.  2  is  the  signal  magnitude  as  a  characteristic 
function  of  contact  impedance  on  exposure  to  radiation  with  X=5  mm. 
As  we  can  see  In  the  figure  presented,  the  polarity  of  the  detected 
signal  varies  by  a  factor  of  2.  In  accordance  with  the  theoretical 
assumptions  elaborated  above,  this  function  can  be  accounted  for  by 
the  effect  of  mirror  image  forces  for  small  film  thicknesses  (small 
contact  impedances)  and  by  the  presence  of  thermionic  emission, 
which  must  be  taken  into  account  for  high  contact  impedances  (large 
film  thicknesses).  When  this  MDM  contact  is  irradiated  with  radia¬ 
tion  with  X*10.6  micrometers,  the  signal  polarity  was  observed  to 
vary  only  for  small  contact  impedances.  Evidently,  the  absence  of 


change  In  polarity  for  large  contact  Impedances  can  be  accounted 
for  by  a  reduction  in  the  effect  of  the  thermionic  emission  current 
when  the  wavelength  is  shortened. 
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Since  the  time  constant  of  the  RC  junction  depends  strongly 
on  the  dielectric  film  thickness,  as  to  be  expected,  for  small 

contact  impedances  the  rise  time  of  the  detected  pulse  front  was 

—7 

small  (less  than  10  s).  As  the  contact  impedance  was  increased, 

-4  -  7 

this  time  rose  to  10  to  10  J  s .  These  functions  were  analogous 
for  both  cases  of  irradiation. 
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PART  4 

dielectric  waveguides 


CALCULATION  OF  MULTILAYER  DIELECTRIC  WAVEGUIDE 

G.  D.  Rozhkov,  A.  S.  Belanov,  and  V.  F.  Vzyatyshev 

Use  of  Dielectric  Waveguides  in  Submillimeter  Wave  Band 

Several  factors  at  present  impede  the  broad  application  of 
dielectric  waveguides  in  the  submillimeter  wave  band.  Of  these, 
the  increase  in  attenuation  and  the  need  to  reduce  the  trans¬ 
verse  dimensions  of  dielectric  waveguides  in  order  to  retain 
the  single  wave  mode  are  leading  barriers. 

Attenuation  in  dielectric  waveguides,  although  slower  than 
in  hollow  metal  waveguides,  intensifies  with  shortening  of  the 
wavelength.  Given  existing  materials,  with  consideration  of 

_ii  , 

losses  of  the  order  of  2*10  ,  for  single  wavelength  waveguides1 

it  reaches  a  value  of  6-10  dB/m,  somewhere  on  the  boundary  between 
the  millimeter  and  the  submillimeter  bands. 

Accordingly,  in  order  to  employ  dielectric  waveguides  in 
the  submillimeter  band,  first  of  all  a  considerable  reduction 
in  attenuation  must  be  attained  in  these  waveguides.  One  possible 
way  of  solving  this  problem  consists  of  using  dielectric  waveguides 
made  of  special  materials.  In  fact,  because  of  features  of  the 
1 

This  refers  to  dielectric  waveguides  with  considerable  delay, 
which  can  be  attained  by  bends  and  twists. 


mechanism  of  propagation  in  dielectric  waveguides  [62],  thermal 
losses  in  the  waveguide  material  is  the  main  cause  of  attenua¬ 
tion,  so  that  attenuation  is  proportional  to  the  material  loss 
angle . 

Is  there  a  possibility  of  reducing  the  loss  angle  of  mater¬ 
ials  in  the  millimeter  and  submillimeter  bands?  Specialists  on 
electrical  properties  of  dielectrics  believe  that  these  possibi¬ 
lities  undoubtedly  exist  [63].  As  early  as  1964,  polyethylene 
with  tan  6=2.7*10-^  at  the  frequency  of  55.2  GHz  [64]  was  pre¬ 
pared.  But  if  this  loss  angle  persists  up  to  the  frequency 
300  GHz,  dielectric  waveguides  made  of  this  material  would  have 
an  attenuation  of  about  1  dB/m  at  a  wavelength  of  1  mm. 

Of  late  there  have  been  reports  about  studies  conducted  by 
several  British  companies  aimed  at  finding  special  materials  for 
dielectric  waveguides  exhibiting  much  lower  losses  in  the  sub¬ 
millimeter  band  than  known  materials.  The  first  achievements 
in  this  direction  are  reported  in  [65]. 

But  even  if  materials  with  loss  losses  were  prepared,  the 
difficulties  Involved  with  the  small  cross-sectional  dimensions 
will  remain.  These  dielectric  waveguides  will  be  difficult  to 
fabricate  with  the  necessary  tolerances  complied  with.  In  addi¬ 
tion,  for  small  dimensions  of  the  dielectric  waveguides  the 
mating  of  waveguide  leads  of  individual  assemblies  when  devices 
are  constructed  of  these  will  be  difficult. 

Specific  difficulties  also  stem  from  the  fact  that  when  any 
external  objects  are  brought  close  to  dielectric  waveguides,  for 
reinforcement,  for  example,  an  external  electromagnetic  field  of 
the  dielectric  waveguide  will  be  excited  and  additional  losses  and 
radiation  will  appear.  Satisfactory  solutions  to  the  problem  of 
reinforcement  have  been  found  in  the  millimeter  band  by  means  of 
members  made  of  foamed  materials,  and  also  with  metal  members. 
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In  the  submillimeter  band,  members  made  of  foamed  materials  can 
become  unsuitable  because  of  the  commensurability  of  their  pore 
size  and  the  wavelength,  while  ohmic  losses  increase  drastically 
in  the  metal  reinforcement  members. 

Since  dielectric  waveguides  as  such,  especially  when  made  of 
polymeric  materials,  exhibit  low  rigidity  and  become  very  diffi¬ 
cult  in  imparting  mechanical  stiffness  to  devices  on  the  dielec¬ 
tric  waveguides,  or  the  stability  of  the  mutual  arrangement  of 
the  dielectric  waveguides  in  sections  with  space _coupling. 

It  is  known  [66]  that  retaining  the  single  wavelength  mode 
requires  that  in  the  dielectric  waveguide  the  so-called  reduced 
dimension 

d  •eH£,-£af*  (l) 

not  exceed  some  value  that  is  critical  for  the  higher  wavelength 
types.  From  Eq.  (1)  it  follows  that  the  physical  dimension  of 
the  dielectric  waveguide  d  can  be  increased  if  the  dielectric 
constant  of  the  dielectric  waveguide  material  is  made  smaller. 
Unfortunately,  the  minimum  value  of  e  of  known  homogeneous  solid 
dielectrics  is  close  to  2.  3ut  using  porous  materials  of  the 
foamed  polystyrene  type  is  hampered  owing  to  the  need  to  produce 
pores  with  a  diameter  much  less  than  the  wavelength,  and  also 
because  the  latter  are  not  rigid  enough. 

However,  there  is  another  possibility  of  increasing  the 
transverse  dimensions  of  the  dielectric  waveguides.  It  consists 
of  increasing  the  dielectric  constant  of  the  external  medium 
in  Eq.  (1).  In  practice,  this  corresponds  to  the  transition 
to  a  multilayer  dielectric  waveguide. 

Multilayer  Dielectric  Waveguides 

A  rod  of  a  material  with  e1  In  a  massive  homogeneous  shell 
of  material  with  e2  is  the  simplest  example  of  this  kind  of 
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of  dielectric  waveguides.  If  this  dielectric  waveguide  is  con¬ 
structed,  for  example,  of  polyethylene  (^=2.28)  in  a  teflon 
jacket  (£2=2.08),  the  transverse  dimensions  of  the  rod  jacket 
can  be  increased  by  roughly  threefold  compared  with  the  case 
of  locating  it  in  air  (e2=l). 

When  a  multilayer  dielectric  waveguide  is  used,  it  is  also 
easy  to  solve  the  reinforcement  problem,  since  the  field  in  the 
jacket  rapidly  fades  with  increasing  distance  from  the  rod.  The 
mating  problem  is  also  eased.  However,  if  the  jacket  surrounds 
the  dielectric  waveguide  core  on  all  sides,  access  to  the  exter¬ 
nal  field  of  the  dielectric  waveguide  is  impeded,  as  is  the 
adjustment  of  the  parameters  of  devices  with  space  coupling. 

We  must  not  forget  that,  in  addition  to  the  useful  effect, 
the  increase  in  the  transverse  dimensions  of  the  guide  rod  and 
the  field  region,  the  directional  properties  are  degraded  in  the 
dielectric  waveguide  jacket.  They  can  be  characterized  by  the 
delay  factor  7i  with  respect  to  the  external  medium,  which  cannot 
be  larger  than  the  quantity  and  in  the  single  wave¬ 

length  mode  generally  it  cannot  exceed  the  quantity  *  ?*,/?  . 

For  example,  in  the  above-discussed  variant  of  a  polyethylene- 
teflon  dielectric  waveguide  }  *  • 


The  poor  directional  properties  of  the  multilayer  dielectric 
waveguide  will  be  exhibited  in  that  radiation  at  the  bent  sections 
will  be  Induced  even  when  the  radii  of  curvature  are  very  large. 

In  this  case,  the  most  significant  fact  may  be  that  the  total 
(for  radiation  and  thermal)  angular  attenuation  for  a  section 
bent  with  the  radius  R 

(2) 

as  shown  in  [66]  cannot  be  smaller  than  some  minimum  quantity 
otg iwm  >  attained  when  there  is  an  optimal  bending  radius  . 
Let  us  evaluate  these  quantities. 
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(3) 


For  angular  attenuation  that  arises  due  to  radiation  in 
dB/rad,  from  [66,673  we  have 

where  J  is  the  delay  of  a  wave  with  respect  to  the  medium  with 
z2 

*'•  %  is  the  relative  radius  of  bending  in  wavelengths  in  the 
ambient  medium  Vvr. 

For  thermal  attenuation  (in  dB/unit  of  length),  from  [66] 
we  can  write 

**  •  W»  ( 4 } 

where  6^  and  6^  are  the  loss  angles  in  the  media  with  and 
k1  and  k 2  are  the  structural  attenuation  coefficients.  In  our 
case ,  ;  here,  k4  '♦  nt • 1  independently  of  the  shape  of 

the  waveguide  cross  section  and  independently  of  its  dimensions. 
If  we  further  set  ,  relation  (4)  becomes 

°*r  *  ^7*^  (b) 

By  substituting  Eqs.  (3)  and  (5)  in  Eq.  (2)  and  differentia¬ 
ting  with  respect  to  R,  for  R‘t  we  get 

C-Wl'i  (6) 

where  y  is  a  root  of  the  equation 

By  substituting  Eqs.  (6)  and  (7)  in  Eq.  (2),  we  have 


(8) 
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Fig.  1 


The  quantities 


(Fig.  a)  and 


(Fig.  b)  as  func¬ 


tions  of  7.  are  presented  in  the  plots  in  Fig.  1.  From  the  plots 

r 

we  see  that,  given  the  existing  materials  (tan  6=2' 10  ),  the 

minimum  losses  in  bending  become  unacceptably  large  even  when 
^<9a?  (larger  than  2  dB/90°).  Here  the  optimal  radius  of  bend¬ 
ing  Is  about  220*2  and  is  obviously  applicable  only  in  the  short¬ 
wave  part  of  the  submillimeter  band  (when  X=0.3  mm  and  £2=2*°^> 
R^-1* 5  mm) . 


Bending  losses  can  be  reduced  considerably  when  tan  5  is 

-7 

made  smaller.  Thus,  when  tan  6=5*10  j  even  the  quantity 

is  applicable,  which  corresponds  to  ftof  and 

makes  It  possible  to  increase  the  transverse  dimensions  by  about 
tenfold  compared  with  the  dielectric  waveguide  in  air.  But  the 
bending  radii  required  in  this  case  (Rl«8-10^)  obviously  are 


-ri| 


unacceptable  through  the  submillimeter  band  (even  when  X=0 .1  mm 
and  e=2.08,  Rmmt  =  5^0  mm). 


Thus,  by  applying  the  simplest  variant  of  a  multilayer  dielec¬ 
tric  waveguide  in  a  situation  when  it  must  be  subjected  to  bends, 
a  substantial  (by  an  order  of  magnitude)  gain  can  be  attained  in 
the  transverse  dimensions  only  in  the  optical  and  in  the  near 
infrared  regions,  and  even  then  under  the  condition  that  the  loss 
angle  of  the  materials  will  be  not  higher  than  10”  to  10~  . 

The  transition  to  the  multilayer  dielectric  waveguide  with 
a  heterogeneous  jacket  opens  another  possibility  of  controlling 
waveguide  parameters.  By  varying  the  dimensions  of  the  indivi¬ 
dual  regions  of  the  jacket  and  the  parameters  of  their  materials, 
we  can  vary  the  properties  of  the  dielectric  waveguide  over  wide 
limits  and  achieve  the  combination  of  materials  that  cannot  be 
achieved  with  a  dielectric  waveguide  with  a  homogeneous  jacket. 

Let  us  refer  to  this  kind  of  multilayer  dielectric  waveguide  with 
a  heterogeneous  jacket  as  a  "composite"  type. 

"Composite"  Dielectric  Waveguides 

A  composite  dielectric  waveguide  is  a  rod  made  of  a  material 
with  ei>E2>£3»  at>uttlnS  on  several  sides  a  material  with  e2,  and 
in  the  remaining  sides,  with  a  medium  that  has  e^,  for  example, 
with  air.  Some  examples  of  these  dielectric  waveguides  are 
satisfactorily  described  with  models  shown  cross-sectionally  in 
Pig.  2. 

As  will  be  shown  below,  the  dimensions  of  the  rod  of  this 
kind  of  dielectric  waveguide,  when  the  single  wavelength  mode  is 
retained,  can  be  nearly  as  large  as  the  rod  dimensions  for  a  homo¬ 
geneous  Jacket  with  e2>  On  the  other  hand,  when  e^»l,  this  kind 
of  dielectric  waveguide  has  a  free  access  to  the  external  field. 

But  the  feature  of  the  behavior  of  the  composite  dielectric 
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waveguide  in  curved  sections  is  the  most  attractive.  Radiation 
from  a  bend  of  this  dielectric  waveguide  must  depend  on  the  direc¬ 
tion  of  the  bend.  In  this  case  it  can  be  expected  that  radiation 
from  a  bend  such  that  the  free  side  of  the  rod  is  the  external 
side  is  determined  by  the  delay  factor  with  respect  to  a  medium 
with  £ ^ .  Then  the  dimensions  of  this  bend  can  be  considerably 
less  than  for  a  dielectric  waveguide  in  a  homogeneous  jacket  with 

e2m 

Unfortunately,  a  rigorous  calculation  of  the  characteristics 
of  the  composite  dielectric  waveguides  is  exceptionally  cumbersome. 
Because  of  additional  interfaces,  this  calculation  must  be  much 
more  complicated  than  the  calculation  of  a  rectangular  dielectric 
waveguide  with  a  homogeneous  jacket.  The  so-called  asymmetric 
plane  dielectric  waveguide,  a  dielectric  layer  with  e.,  a  thickness 
2a,  and  whose -medium  on  one  side  has  the  dielectric  constant 
and  on  the  other,  e2,  is  a  relatively  simple  model  satisfactorily 
describing  this  dielectric  waveguide  when  4»*o  . 


Analysis  of  the  asymmetric  plane  dielectric  waveguide  is  re¬ 
ported  by  the  studies  [68,69].  We  carried  out  a  numerical  solu¬ 
tion  of  the  characteristic  equation  and  calculated  the  wave  numbers, 
the  constant  of  propagation,  the  effective  field  dimensions  on  both 
sides  of  the  wafer,  and  the  structural  attenuation  coefficients  due 
to  losses  in  the  wafer  material  and  in  the  adjacent  media  for  the 
waves  H1,  H2,  H^,  and  E2,  E^,  and  the  following  parametric 
ratios : 
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;  f, -2,28;  £,  -I;  £,-1;  I.'oS;  1.3;  and  2 . 08 ; 

*>3,9;  £>  I;  £a- I;  1.05  and  2.08  over  a 
wide  range  of  values  of  the  reduced  wafer  thickness. 

Lacking  the  possibility  of  presenting  all  the  results  in 
this  present  study,  we  give  just  the  dispersion  characteristics 
for  the  case  ^=2.28  (polyethylene);  £^=1;  a2=2*0^  (teflon). 

They  are  noted  with  the  symbol  ••in  the  plots  in  Pigs.  3  and  4. 

The  upper  curves  correspond  to  the  E^  wave  in  Pig.  a  and  to  the 
wave  in  Pig.  b,  and  the  lower  curves,  to  the  E2  and  H2  waves. 
The  quantity  U  in  the  plots,  the  delay  factor  of  the  wave  with 
respect  to  the  medium  with  as  associated  with  the  earlier 
discussed  quantity  by  the  relationship 

(9) 

As  can  be  seen  from  the  plots,  the  quantity  U  satisfies  the  in¬ 
equalities 

(10) 

here  executing  the  leftmost  equality  in  Eq.  (10)  corresponds  to 
the  critical  frequency,  and  of  the  rightmost  equality,  to  infi¬ 
nitely  large  frequencies.  The  physical  meaning  of  the  leftmost 
inequality  in  Eq.  (10)  corresponds  to  the  case  when  the  wave  field 
in  a  medium  with  e2  decays  exponentially  with  increasing  distance 
from  the  guide  rod  with  e  ;  when  it  is  not  satisfied,  the  field  in 
a  medium  with  e2  will  taken  on  the  form  of  plane  waves  propagating 
at  some  angle  to  the  z-axis  and  entraining  energy  from  the  guide 
rod. 
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a*  a  to  U 

b) 

Pig.  3 
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Pig.  4 


! 
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Relative  Dimensions  and  Band  Dimensions  of  an  Asymmetric 
Plane  Dielectric  Waveguide 

Analysis  of  the  critical  conditions  for  the  nearest  highest 
wave  type  shows  that  the  wafer  thickness,  at  parameters  character¬ 
istics  of  the  curves  in  Figs.  3  and  4,  can  be  chosen  at  values 
that  are  larger  by  a  factor  of  3.5  to  3.6  than  the  thickness  of 
the  wafer  in  air  (^=2.28;  e2=l).  Interestingly,  in  the  case  of 
a  homogeneous  setting  (£^*2.28;  £2=£^=2.08)  an  increase  in  wafer 
thickness  of  only  a  factor  of  2.5  (by  30  percent  less)  is  possible. 

It  is  interesting  also  to  compare  the  band  properties  of  these 
dielectric  waveguides.  Considering  that  the  minimum  wavelength  of 
the  band  A,**,  corresponds  to  the  critical  conditions  of  the  nearest 
highest  wave  type  and  assuming  that  at  the  maximum  wavelength,  the 
delay  must  be  equal  to  some  minimum  value  ,  for  the  band 

overlap  factor  we  have 


Assuming  in  this  case  J>Mmm  =0 . 02 ,  for  the  asymmetric  dielectric 

waveguide  we  have  D„  *1.67  and  D-,  *1.58,  while  for  the  symmetric 

nl  ^1 

variant  Du  *1.56  and  Dt,  *1.49  (9  to  11  percent  less). 

H1  E1 

Thus  the  plane  asymmetric  dielectric  waveguide  will  exceed 
the  symmetric  variant  both  with  respect  to  the  allowable  trans¬ 
verse  dimensions  and  with  respect  to  the  width  of  the  working  band. 

Possibilities  of  Approximate  Calculation 

To  calculate  the  models  In  Fig.  2  that  are  closer  to  the  actual 
situation,  we  can  employ  an  approximate  approach  first  described  by 
V.  Shlosser  in  [70]  and  widely  used  by  Ye.  Markatili  In  [73].  It 
is  described  in  detail  also  in  the  review  [72].  This  approach  is 
based  on  the  assumption  that  the  fields  in  the  regions  /x(«a  and 
lyl  * S  can  be  described  with  a  single  trigonometric  function  within 


the  rod  and  by  a  single  exponential  function  outside  the  rod,  and 
with  fields  in  the  region  /x/»d;  /y/*>  6  need  not  generally  interest 

us.  In  [71]  it  is  shown  that  this  approach  for  small  e.  yields 
satisfactory  results,  especially  in  the  region  that  is  far-removed 
from  the  critical  frequency. 

The  characteristic  equation  derived  by  the  method  described 
can  be  represented  in  the  form 


where  0  is  a  dimensionless  internal  wave  number  [66]  of  a  plane 
dielectric  waveguide  with  thickness  2b  made  of  a  material  with 
in  a  medium  with  £2*  the  case  of  model  1,  and  in  a  medium 
with  e 2 ,  in  the  case  of  model  2;  and  is  the  delay  coefficient 
(with  respect  to  a  medium  with  £^)  of  the  above-examined  plane 
asymmetric  dielectric  waveguide. 

The  results  of  calculating  the  delay  with  respect  to  Eq.  (12) 
for  dielectric  waveguides  with  different  b/a  ratios  are  presented 
in  Pigs.  2  and  4. 

The  symbols  for  both  wave  types,  as  In  [72],  are  as  follows: 
the  letter  symbol  [E]  gives  the  field  (in  this  case,  electric), 
with  respect  to  which  the  polarization  direction  is  estimated; 
the  letter  superscript  (x  or  y)  gives  the  coordinate  axis  in 
whose  direction  the  indiicated  field  is  generally  polarized;  the 
numerical  subscripts  give  the  number  of  field  variations  in  the 
directions  of  the  corresponding  coordinate  axes. 

It  is  to  be  expected  that  the  accuracy  of  these  results  im¬ 
proves  as  U  grows  larger.  The  values  (dashed  sections  of 

the  curves)  will  contradict  the  physical  considerations,  since 
the  waves  traveling  with  a  delay  v<irt  must  radiate. 

Estimates  of  the  width  of  the  working  band  based  on  criterion 
(11)  show  that  it  decreases  with  decrease  in  b/a.  When  b/a*2,  the 


decrease  is  not  large  and  amounts  to  2-7%.  When  b/a=0.5»  model  2 
is  a  much  more  advantageous  model  from  the  standpoint  of  band 
properties . 


EXPERIMENTAL  STUDY  OF  A  COMPOSITE  DIELECTRIC  WAVEGUIDE 
V.  F.  Vzyatyshev,  G.  D.  Rozhkov,  and  B.  A.  Ryabov 
Formulation  of  Problem 

When  any  microwave  transmission  line  is  investigated,  it  is 
necessary  first  to  establish  the  conditions  for  single-wavelength 
or  a  virtually  single-wavelength  operating  conditions  and  to  study 
the  main  characteristics  (delay,  attenuation,  and  distribution  of 
field)  of  the  waves  under  these  conditions.  A  method  for  the  ex¬ 
perimental  resolution  of  these  problems  as  applied  to  dielectric 
waveguides  has  been  worked  out  (see  C73>7^]). 

One  feature  of  setting  up  experiments  with  a  composite  dielec¬ 
tric  waveguide  consists  of  the  fact  that  at  the  present  time  there 
are  still  numerous  assemblies  (exciters,  bends,  and  so  on)  that  are 
needed  for  known  methods  of  measuring.  So  we  were  compelled  to 
work  out  a  special  method  of  measuring,  using  the  minimum  number 
of  assemblies  that  are  simplest  in  a  design  sense. 

At  the  first  stage  we  set  up  two  main  goals. 

1.  Investigation  of  the  dispersion  characteristics  of  the 
lower  wave  types. 


2.  Investigation  of  the  behavior  of  these  waves  at  the  bend 
sections . 


Fig.  1 


The  cross  section  of  the  waveguide  investigated  is  shown  in 

Fig.  1.  Rod  1  is  made  of  polyethylene  (e  =2.28)  and  has  the  dimen- 
p  -*■ 

sions  10x4.5  mm  .  It  is  a  rectangular  PDV-12  dielectric  waveguide 
that  is  truncated  to  make  one  of  its  faces  planar  [73].  Region  2 

p 

is  made  of  teflon  ( £2=2.08)  and  has  the  dimensions  10x40  mm  . 

This  waveguide  configuration  is  advantageous  in  that  it  provides 
for  access  to  the  external  field  of  the  dielectric  waveguide  from 
three  sides. 

Dispersion  Characteristics 

This  investigation  was  conducted  with  a  dielectric  resonator 
and  following  the  method  described  in  [73 >7*0,  in  the  wavelength 
band  from  4  to  10  mm.  Serving  as  the  resonator  was  a  section  of 
the  above-described  composite  dielectric  waveguide,  120  mm  long, 
placed  between  a  solid  and  a  gridded  reflector.  The  resonator 
is  excited  through  the  gridded  reflector  with  the  open  end  of  a 
rectangular  dielectric  PDV-8  waveguide,  with  a  cross  section  of 
6.7x3.35  mm  .  The  back  wave  is  branched  with  a  balanced  wave¬ 
guide  joint  on  this  same  waveguide. 

It  was  found  that  when  the  exciting  dielectric  waveguide  is 
placed  opposite  the  midsection  of  the  teflon  plate,  a  large  number 
of  resonances,  partially  superimposing  on  each  other,  is  observed. 


jjuiuw.jww-y 


Identification  of  the  type  of  these  oscillations  was  made 
difficult  owing  to  the  high  concentration  of  their  field  in  the 
dielectric.  When  the  exciting  dielectric  waveguide  is  placed 
opposite  the  polyethylene  rod,  the  intensity  of  most  resonances 
is  reduced,  and  the  intensity  of  resonances  corresponding  to  the 
wave  types  we  are  interested  in  is  even  increased. 


The  longitudinal  order  of  the  resonance  oscillation,  N,  the 
number  of  half-waves  accommodated  along  the  resonator  length  L, 
is  determined  with  an  electrical  probe,  and  the  resonance  frequency 
f  ,  with  a  wavemeter.  Delay  with  respect  to  air  is  then  calculated 
based  on  the  relationship 


(1) 


where  C  is  the  speed  of  light. 


Polarization  of  the  electric  field  of  the  observed  oscilla¬ 
tions  and  the  number  of  Its  variations  are  determined  with  a  thin 
plate  made  of  metallized  mica. 


The  experimental  results  are  shown  with  curves  marked  by 
circles  in  Fig.  2.  Also  shown  there  are  the  calculated  results 
obtained  with  the  approximate  method  described  in  [118]  (solid 
curves  in  the  upper  part  of  the  plot).  The  symbols  for  the  wave 
types  correspond  to  those  in  [118].  Clearly,  the  discrepancy 
of  the  experimental  results  and  the  calculation  Is  very  substan¬ 
tial.  Admittedly,  the  experimental  waveguide  (Fig.  1)  differs 
from  the  calculation  models  (Fig.  1  in  [1183) ,  but  this  difference 
Is  generally  not  taken  into  account  in  the  approximate  calculation 
method. 


The  most  characteristic  experimental  fact  is  the  detection  of 
waves  with  .the  delay  This  fact,  which  at  first  glance  con¬ 

tradicts  the  physical  principles  advanced  in  [118],  is  easily  ex¬ 
plained  if  we  consider  that  the  region  with  e2  has  a  finite  height 
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Pig.  2 


(10  ram).  So  the  waves  entraining  energy  from  the  polyethylene  rod 
at  frequencies  lower  than  the  critical  frequency  must  be  not  plane 
waves  in  an  unbounded  medium,  but  waves  in  the  dielectric  plate 
(e2=2.08;  2b=10  mm).  The  dispersion  characteristics  of  these 
waves  are  indicated  with  solid  curves  in  the  lower  part  of  the 
plot  in  Fig.  2.  The  symbols  for  these  waves  correspond  to  those 
in  [66].  From  the  plots  we  see  that  some  experimental  points  are 
below  these  curves,  although  from  physical  considerations  it  follows 
that  they  always  must  be  higher  and  can  be  approached  only  in  the 
critical-frequency  region.  But  this  fact  cannot  be  accepted  as 
proven,  since  the  discrepance  falls  within  the  limits  of  experi¬ 
mental  error  (about  2 *10”^),  the  calculation  (about  5*10  ^),  and 
the  dielectric  constant  of  teflon  (about  5*10  ^). 

Behavior  of  Waves  at  Bend  Sections 

The  measurements  were  made  in  the  wavelength  band  of  3-9  to 
5.7  nun  using  the  following  method. 


Pig.  3a 


The  transmission  coefficient  of  the  waveguide  line  consisting 
of  a  1  m  long  section  of  a  polyethylene  PDV,  with  transverse  dimen¬ 
sions  of  10x4.5  mm,  the  midsection  of  which  can  be  converted,  by 
bringing  up  a  rectangular  teflon  plate  (Fig.  3a)  or  shaped  teflon 
plates  (Fig.  3b),  into  a  section  of  a  linear  or  a  bent  composite 
dielectric  waveguide  and  two  exciters  in  the  form  of  smooth  wave- 
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guide  metal  tapers  with  a  3. 6x1. 8  mm  cross  section  for  a  section 

2 

with  the  dimensions  of  5.5x11  mm  and  a  100  mm  length.  The  entrance 
and  exit  sections  of  the  taper  are  turned  by  90°  in  order  for  the 
electric  field  in  the  output  section  to  be  oriented  parallel  to  its 
large  side. 


Losses  in  a  line  made  up  of  a  PDV  with  tapers  for  a  wavelength 
of  3.9  mm  came  to  about  6  dB.  With  the  thermal  losses  in  the  dielec¬ 
tric  waveguide  material  and  in  the  tapers  subtracted,  the  losses  in 
excitation  are  about  1.2  dB  per  exciter.  As  the  wavelength  is  made 
larger,  the  losses  are  reduced.  In  order  to  reduce  the  relative 
amplitude  of  the  higher-type  waves,  special  filters  are  inserted 
near  the  exciters.  By  applying  an  attenuation  of  about  5dB  on  the 
principal  wave  type,  the  filters  serve  simultaneously  for  decoupling 
the  members  being  measured,  the  generator  and  the  detector. 

It  was  established  that  the  PDV  bends  in  the  plane  of  its 
smaller  facet  with  radii  to  6  cm  for  an  angle  to  l8o°  do  not  intro¬ 
duce  additional  losses  of  more  than  0.8  dB.  The  losses  observed 
do  not  depend  monotonically  on  the  angle  of  bending.  Probably  they 
are  associated  with  the  conversion  to  higher-type  waveguide  waves. 

The  main  result  of  the  measurements  is  the  losses  when  shaped 
teflon  plates  (Figs.  3b  and  c)  are  brought  up  to  the  bent  section 
of  the  PDV.  This  section  is  thus  converted  to  a  bent  section  of 
a  composite  dielectric  waveguide  (a  bend  of  the  kind  shown  in 
Fig.  3b,  when  the  material  with  e2  is  in  the  external  side  of  the 
bend,  we  call  an  external  bend,  and  a  bend  of  the  Fig.  3c  type,  an 
Internal  bend).  It  was  found  that  the  losses,  especially  in  the 
external  bend,  depend  strongly  on  the  gap  between  the  PDV  and  the 
plate.  So  measures  were  adopted  to  reduce  the  gap. 

Since  during  transit  through  the  section  with  the  plate  at 
the  mating  of  the  dielectric  waveguide  with  the  composite  dielec¬ 
tric  waveguide  there  must  be  some  conversion  into  both  the  higher- 
type  waves  and  into  radiation  waves,  the  characteristics  of  this 
converter  were  measured  separately.  To  do  this,  the  transmission 
coefficient  was  plotted  as  a  function  of  the  distance  between  the 
PDV  and  the  Fig.  3a  plate  brought  up  to  the  PDV,  so  that  on  touch¬ 
ing,  the  composite  dielectric  waveguide  of  Fig.  1  was  constructed. 
This  function  in  all  cases  proved  to  be  monotonic,  with  the  addi¬ 
tional  losses  not  exceeding  0.5  dB,  X«5.7  mm,  and  0.2  dB,  A*3.9  mm. 
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Additional  losses  in  the  internal  bend  cannot  be  measured 
exactly  because  during  the  bringing  up  and  pressing  up  of  shaped 
plates,  the  position  of  the  PDV  ends  in  the  exciters  and  the  losses 
in  excitation  are  changed.  However,  it  can  be  ensured  that  addi¬ 
tional  losses  do  not  exceed  1  dB  at  the  wavelength  of  5.7  nun  and 
even  less  so  at  shorter  wavelengths. 

In  the  external  bend,  as  to  be  expected  [118],  the  losses 
prove  to  be  much  larger.  They  increase  with  increase  in  wavelength 
and  amount  to  20  dB  at  1=3.9  mm  and  to  35  dB  at  1=5.7  mm. 

Thus,  the  composite  dielectric  waveguide  with  heterogeneous 
surrounding  media  actually  has  certain  advantages  when  compared 
with  the  multilayer  dielectric  waveguide  In  a  homogeneous  setting. 


DIELECTRIC  WAVEGUIDE-BASED  WAVE  TYPE  CONVERTER 

A.  N.  Merkur'yev 

Approach  to  Converter  Analysis 

As  we  know,  when  another  single  wavelength  dielectric  waveguide 
brought  up  to  the  wave  field  propagating  along  a  first  single  wave¬ 
length  dielectric  waveguide,  energy  branches  from  the  first  to  the 
second.  Phenomena  in  this  system  can  be  described,  with  an  accu¬ 
racy  suitable  for  a  number  of  practical  cases,  with  the  theory  of 
bound  waves  [75].  It  is  shown  in  [76]  that  the  level  of  energy 
interchange  between  the  dielectric  waveguides  is  determined  by  the 
ratio  2c/Ay,  where  c  is  the  linear  coupling  coefficient  (LC)  between 
the  waves,  and  Ay  is  the  difference  of  their  constants  of  propaga¬ 
tion.  If  Ay=0,  there  may  be  total  branching  of  energy  from  the 
initially  excited  dielectric  waveguide  into  the  second  waveguide; 
this  will  occur  when 

cL  •  f  *  **,  (  n  •  <J  ( l ) 

where  L  is  the  coupling  section  length. 

But  if  Ay?<0  and  2c/Ay<<l,  as  follows  from  [76,79],  the  frac¬ 
tion  of  energy  branching  into  the  second  dielectric  waveguide  will 
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not  exceed  (2c/Ay)  .  This  fact,  first  pointed  out  in  [75]  as  applied 


to  hollow  metal  waveguides,  then  can  be  used  in  converting  the 
principal  wave  of  the  single  wavelength  dielectric  waveguide, 
for  example,  a  multiwavelength  dielectric  waveguide,  in  particular, 
to  an  HQ1  wave  of  the  circular  dielectric  waveguide. 

Analysis  of  phenomena  in  the  system  of  coupled  single 
wavelength  and  multiwavelength  dielectric  waveguides  is  hampered 
by  the  fact  that,  since  a  single  wavelength  transmission  line  must 
be  brought  into  correspondence  with  each  propagating  wave  type, 
the  number  of  differential  equations  in  their  system  describing 
the  coupling  of  the  dielectric  waveguides  is  relatively  large, 
so  that  numerical  methods  must  b'e  used  in  the  solution.  The  ana¬ 
lysis  can  be  simplied  with  allowance  for  the  fact  that  the  wave 
type  converter  must  operate  with  minimum  losses,  that  is,  in  it 
there  must  be  a  substantial  coupling  of  the  wave  from  the  single 
wavelength  dielectric  waveguide  only  with  a  single  wave  type  from 
the  multiwavelength  dielectric  waveguide;  coupling  with  other  wave 
types,  as  an  unwanted  phenomenon  degrading  the  device  characteris¬ 
tics,  must  be  brought  to  a  minimum. 

Obviously,  this  can  be  achieved  by  selecting  the  cross-sec¬ 
tional  dimensions  and  the  dielectric  constant  of  the  dielectric 
waveguide  so  that  there  is  an  equality  of  the  constants  of  propa¬ 
gation  of  the  converted  wave  types,  and  by  reducing  the  linear 
coupling  coefficient;  this  requires  making  the  distance  between 
the  dielectric  waveguides  quite  large  [771. 

Further,  to  describe  the  phenomena  in  the  dielectric  wave¬ 
guide  system  we  can  use  the  approach  presented,  for  example,  in 
[78],  that  is,  we  can  simultaneously  examine  the  coupling  of  only 
two  specified  wave  types  without  allowing  for  the  transformation 
into  parasitic  wave  types,  and  the  latter  can  be  estimated  by  the 
weak  coupling  method.  To  do  this,  the  results  obtained  in  [79]  can 
be  used  directly  in  calculating  the  converter  characteristics; 


this  however  presupposes  a  knowledge  of  the  linear  coupling  coeffi¬ 
cient.  Since  it  is  not  possible  to  calculate  the  linear  coupling 
coefficients,  we  calculate  them,  more  precisely,  the  linear  coupling 

coefficients  of  the  principal  HE,,  wave  of  the  PDV  and  the  H  wave 
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oi  the  XDV,  experimentally. 

Measuring  the  Linear  Coupling  Coefficients  of  Waves  HE^V  and 

The  experimental  estimate  of  the  linear  coupling  coefficients 

is  conveniently  done  by  investigating  the  characteristics  of  an 

adjustable  converter  mockup.  The  measurements  were  made  according 

to  the  diagram  in  Pig.  1.  The  converter  I  as  such  consists  of  a 

XDV  section  4.35  mm  in  diameter  and  a  PDV  section,  with  a  cross 
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section  of  2.90x1.45  mm  ,  parallel  for  a  coupling  section  280  mm 
long.  These  measurements  of  the  dielectric  waveguide  sections 
provide,  as  is  clear  from  the  dispersion  curves  in  Fig.  4,  an  equa¬ 
lity  of  delays  £=A/Ah-1  (A.  is  the  wavelength  in  the  dielectric 

D  ^  PDV  KDV 

waveguide)  for  the  waves  HEn  anc*  Hoi  in  the  wavelenSth  band 
A%5.1  mm.  In  this  case  (see  Fig.  2),  in  addition  to  the  principal 
wave  HE^-j-  and  the  wave  HQT,  two  more  wave  types  can  be  propagated 
in  the  XDV:  EQ1  and  HE21«  Coupling  with  the  first  of  these,  EQ1, 
is  eliminated  by  selecting  the  mutual  positioning  of  the  dielec¬ 
tric  waveguides  such  as  in  shown  in  the  field  in  Fig.  3*  In  the 
Fig.  1  diagram  the  reflection  coefficient  of  the  converter  is 
measured  (reflector  2  is  connected  to  the  XDV  of  this  converter). 

The  following  adjustments  are  possibled  in  the  converter:  varying 
the  distance  t  between  the  dielectric  waveguides  (see  Fig.  3)  and 
shifting  one  dielectric  waveguide  along  the  other  with  constant  t. 
Both  adjustments  were  made  without  changing  the  length  of  the  coup¬ 
ling  section. 

As  shown  by  the  experiment,  when  one  dielectric  waveguide  is 
shifted  along  the  other,  the  reflection  coefficient  rapidly  oscil¬ 
lates;  a  slow  aperiodic  change  is  superimposed  on  the  oscillation. 

The  latter  is  evidently  caused  by  some  disturbance  in  the  parallelity 
of  the  dielectric  waveguides  when  they  are  shifted.  Rapid 
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oscillations  are  explained  by  the  interference  in  the  detector  of 
two  (or  possibly  several)  reflected  waves.  It  turned  out  that  near 
1=5.1  mm,  when  t>2.5  mm,  the  period  of  oscillations  corresponds  to 
the  delay  of  the  wave  in  the  KDV.  Therefore,  in  the  converter 
there  must  actually  be  the  predominant  coupling  with  the  wave 
in  the  KDV,  and  the  oscillations  are  caused  by  the  interference 
of  the  wave  with  the  variable  phase  reflected  from  the  converter, 
and  the  wave  with  a  constant  phase,  caused  by  a  nonideal  matching 
of  the  arrangement  elements. 


Pig.  1 

1.  adjustable  converter  6.  attenuator 

2.  reflector  7.  generator 

3,4.  tapers  to  the  PDV  from  8.  recording  microammeter 

metal  waveguides  9.  amplifier 

5.  bridge  (3  dB)  for  the  PDV 


The  reflection  coefficient  as  functions  of  the  distance  t 
between  the  dielectric  waveguides  in  accordance  with  [77]  is 
oscillatory.  Applying  the  method  given  in  [77],  from  these  func¬ 
tions  we  found  the  linear  coupling  coefficients  of  the  waves 
HEi?V  and  H01V>  whose  Plots  are  shown  in  Fig.  3. 


1.  KDV 

2.  PDV 

3.  supports 


Pig.  2 


4.  housing 

5.  tuning  screw 


Construction  and  Characteristics  of  Converter 

Further  study  of  a  converter  of  HE^V  and  Hq^V  wave  types 
was  made  with  a  nonadjustable  mockup  (see  Fig.  2).  In  it  use 
was  made  of  the  dielectric  waveguides  of  the  same  transverse 
dimensions  and  the  same  waveguide  positioning  as  in  the  adjustable 
converter  mockup.  The  dielectric  waveguides  were  mounted  in  foam 
plastic  supports.  The  ends  of  the  dielectric  waveguides  were 
smoothly  diverted  from  each  other. 


In  trying  to  attain  small  device  dimensions,  we  selected  the 
distance  between  dielectric  waveguides  t  =  2  mm.  Further,  on  the 
basis  of  Fig.  3  and  condition  (1)  for  M=0,  we  found  the  length  of 
the  coupling  section,  L=90  mm. 

The  shape  of  the  PDV  axis  in  the  coupling  section  can  be 
modified  somewhat  with  a  tuning  screw,  attaining  the  maximum  trans¬ 
mission  coefficient. 

In  tuning,  the  two  converters  were  connected  in  a  cascade 
manner.  By  rotating  the  tuning  screws,  the  minimum  losses  A 
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was  attained,  defined  as  the  ratio  of  the  input  signal  strength  to 
the  output  signal  strength.  The  plot  of  the  minimum  losses  A  min 
for  a  single  converter  as  a  function  of  wavelength  is  presented 
in  Pig.  5.  The  quantity  ^min  has  its  smallest  value  at  X=4.7  mm, 
and  not  at  X=5.1  mm,  as  would  be  expected  on  the  basis  of  Pig.  4. 
The  concern  evidently  is  that,  for  total  interchange  of  energy,  it 
is  required  to  have  an  equality  of  velocities  of  propagation  of 
the  coupled  waves;  this  occurs  at  the  point  of  the  frequency  scale 
not  coinciding  with  the  equality  point  of  the  phase  velocities  of 
the  waves  in  the  connected  dielectric  waveguides;  and  it  does  so 
the  more  strongly,  the  closer  the  coupled  dielectric  waveguides 
are  to  each  other. 
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Pig.  5 

When  the  frequency  characteristic  was  being  recorded,  as  before 
the  two  converters  were  connected  cascadewise.  The  losses  for  a 
single  converter  as  a  function  of  wavelength  A(X)  during  the  tuning 
of  the  converters  at  XQ=4.7  mm  is  shown  in  Pig.  6.  In  the  region 
X=»4.25  mm  we  observed  a  second  hump  of  the  frequency  characteristic; 
this  apparently  is  caused  by  the  amplification  of  the  coupling  of 
the  PDV  wave  HE^  and  of  the  KDV  wave  ^21’  The  concern  is  that, 
as  shown  experimentally,  even  at  the  wavelength  4.7  mm,  the  principal 


PDV  wave  is  partially  converted  to  an  HE^  wave  existing  in  the 
KDV  together  with  the  Hni  wave,  and  with  a  reduction  in  X  ,  the 
weakening  of  the  coupling  of  the  HE^  wave  with  the  wave 

is  partially  compensated  by  an  amplification  of  the  coupling  of 


the  HE^V  and  the  waves. 


Thus,  the  minimum,  losses  occur  at  the  wavelength  of  4.7  mm 
and  amount  to  2  dB.  In  the  wavelength  band  from  4.6  mm  to  4.85  mm 
the  losses  do  not  exceed  3  dB.  In  this  mockup,  conversion  of  the 
KDV  wave  to  the  parasitic  type  was  noted. 


Pig.  6 
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